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ABSTRACT 

Solar-cell art ays provide the primary energy source for most 
earth-orbiting satellites. The solar cell utilizes the radiative en- 
ergy from the sun in direct photovoltaic conversion to provide elec- 
trical power to the spacecraft. The computer program described 
here will account for all factors - including isolation diode voltage 
losses, charged particle irradiation, temperature, solar incidence 
angle and other solar cell degradation factors - affecting the power 
output of solar arrays, and thereby greatly reduce the lengthy and 
laborious hand calculations now necessary In designing an array or 
predicting its performance. 

The program will determine the total power output oi solar ar- 
rays consisting of up to 25 electrically-parallel panels; each panel 
may have its own series-parallel solar cell arrangement, temper- 
ature verses time profile and sun angle verses time profile. The 
output of the program is a series of current-voltage (I-V) points 
representing the total solar array I-V characteristic at up to 20 
times during an orbit. 
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THE SOLAR ARRAY SYNTHESIS COMPUTER PROGRAM 


I. INTRODUCTION 

The accurate determination of a solar array's power output in orbit is es- 
sential in the adequate design of a spacecraft power system. The lengthy and 
laborious hand calculations now associated with this determination can be greatly 
reduced through the use of the Solar Array Synthesis Computer Program. This 
program will compute and print out solar array output I-V (current-voltage) 
characteristics for arrays consisting of up to twenty-five (25) solar cell panels. 
Each panel may have its own series - parallel solar cell arrangement, temper- 
ature verses time and solar incidence angle verses time profiles. Also, the 
program will account for isolation diode voltage losses , charged particle ir- 
radiation degradation and other solar cell degradation factors. Calculations 
may be performed for cells of either 1 or 10 ohm-cm nominal base resistivity. 

An individual solar cell I-V characteristic is entered into the computer 
data deck as a series of current and voltage points at nominal solar intensity 
and room temperature . After calculating the equivalent 1 -Mev electron flux from 
the given charged particle spectrum, the program will degrade the solar cell 
characteristic to account for the equivalent 1-Mev flux. The cell I-V curve is 
next shifted to account for user-specified current and voltage degradation factors. 
The fully-degraded ,-solar cell I-V characteristic is then extrapolated into a 
family of solar cell I-V curves at nominal intensity for 15 temperatures (covering 
the temperature range of interest to the user) through the use of specified cur- 
rent and voltage temperature coefficients and appropriate changes to the I-V 
curve shape. Each panel’s contribution to the total solar array I-V curve is de- 
termined separately depending upon the panel’s solar cell arrangement, temper- 
ature and sun angle; all solar cells on any one panel must be operating at the same 
temperature and solar incidence angle. 

Summing the individual panel output currents at a voltage as voltage is 
increased from zero to the open-circuit value of the highest-voltage panel, in 
one-volt increments, yields the total solar array I-V characteristics. Up to 20 
such solar array I-V curves can be obtained during a single program execution, 
corresponding to 20 equal time increments during an orbit, or 20 combinations 
of solar panel temperature and illumination incidence angle as defined by the 
program user in the input data. 

The program consists of five parts: 

• MAIN - reads input data, prints output data, calls the subroutines as 
required, computes and stores individual panel I-V curves for up to 


twenty-five panels and computes the total resulting solar array I-V 
curve in one volt increments . 

• PHI - computes a damage - equivalent, normally - incident (deni) 1- 
Mev electron flux from tables of omnidirectional electron, proton, 
solar flare proton and solar flare alpha particle fluxes , accounting for 
the protection afforded by the selected values of coverglass and back- 
shielding thicknesses. 

9 DEGRAD - degrades a room temperature solar cell I-V curve for the 
equivalent 1-Mev electron flux computed by Subroutine PHI or any other 
user-specified value of equivalent 1-Mev electron flux. 

• STASH - degrades a solar cell I-V curve for current and voltage degra- 
dation factors and series resistance effects and computes and stores 
the degraded cell I-V curve over the specified temperature range. 

• STINT - stores in tables, and then supplies on command, the values of 
variables which are functions of one, two or three arguments. The name 
stands for Standard Table iNT erpolation, and is an adaption of the STINT 
routine which is a part of the IBM SHARE library. 

The program user may elect to perform solar cell calculations only, or 
solar cell and solar array calculations . Any number of runs may be chained 
together if desired. The program, originally written in the FORTRAN II com- 
puter language compatible with the IBM 7094 computer, has been converted to 
the FORTRAN IV language for use on the IBM 360 computer. The information 
contained in this report pertains equally to both versions. 


H. PROGRAM DESCRIPTION 

This section provides a description of each of the program subroutines . 
Figure 1 presents a basic block diagram of the solar array synthesis program 
and summarizes the important functions of each subroutine. 

A. MAIN - Computer Program Control 

The primary function of the MAIN routine is to perform all data input and 
parameter initialization, and call the other subroutines as required. In addition, 
MAIN performs all of the solar array calculations . Up to 20 solar array I-V 
curves can be determined during one computer run, corresponding to 20 equal 
time increments during the sunlight-illuminated portion of an orbit or any 20 
specified combinations of solar panel sun angle and temperature. 
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Figure 1 . Subroutine Functions in Solar Array Synthesis Computer Program 
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The contribution to the total solar array I-V curve of each of up to 25 
parallel panels is calculated Individually, After first determining the panel’s 
operating temperature and solar incidence angle, a representative solar cell 
curve is constructed, A total panel I-V curve is determined by multipling the 
representative cell’s current values by the number of parallel strings on the 
panel, and the voltage points by the panel's number of series solar cells. (Each 
parallel string of a given panel must have the same number of series cells, 
although separate panels may havo different series - parallel combinations . ) If 
the program user has specified a value of array blocking diode loss, this drop 
is incorporated into the calculation of the panel I-V curve. Summation of the 
individual panel curves yields a total solar array I-V characteristic at each timo 
increment. 

D. Subroutine PHI - Calculation of 1-Mov Electron Flux 


Subroutine PHI converts the known charged particle spectrum of the four 
types of particles which damage solar colls most: trapped orbital electrons, 
trapped orbital protons, solar-flare protons and solar-flare alpha particles, to 
a damage-equivalent, normally-incident (deni) 1-Mov electron flux. Solar cell 
irradiation degradation is usually defined in terms of this low-energy flux be- 
cause the damage equivalency of 1-Mev electrons for each type and energy range 
of space particle has been defined. 1 Once calculated, the 1-Mev equivalent flux 
for each type of particle is summed to determine the total 1-Mev flux the cell 
experiences . 

A flux table, prepared for each type of particle, contains the particle’s 
flux population through specific energy ranges for a particular time in orbit. 
Tables 1 through 4 list the energy ranges for each type of particle and gives the 
population of each at one year in a given example orbit. The electron environ- 
ment is broken into 11 ranges; the other particle environments are divided into 
18 ranges. The program will operate properly only when the particle populations 
are divided into the energy ranges specified in tables 1 through 4 . 

The program calculates the equivalent 1-Mev electron flux for each type 
of particle by applying a damage factor (KD) to the population of each energy 
range; damage factors depend on the effective coverglass and backshielding 
densities (gm/cm 2 ) protecting the solar cell. Tables 5 and 6 list damage factors 
for orbital electrons and protons, based on the data reported by Brown, Gabbe 
and Rosenzweig. 2 Although previously trapped orbital protons were distinguished 
from solar-flare protons-based primarily on the relative uncertainty associated 
with the solar-flare proton model-both proton classifications will use the same 
damage factor table . Damage factors for alpha particles , from the work of 
Smith and Blue, 3 are based on the relationship: 

3.8* p (E) = * a (4E) 
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Tatyle 1 

Electron Flux Format 


Energy Range Number 



AE (Mev) 


0.0-0.25 
0.25-0.50 
0.50-0.75 
0.75-1.0 
1 . 0 - 2.0 

2. 0- 3.0 

3. 0- 4.0 
■4. 0-5.0 

5. 0- 6.0 

6. 0- 7.0 
>7.0 


Ne (Electrons/cm 2 ) 


1.06E15 

8.58E14 

2.84E14 

9.12E13 

3.87E12 


1.49E11 

1.42E9 

1.08E7 

1.13E5 

6.57E4 

3.65E3 


Table 2 

Proton Flux Format 


Energy Range Number 



A E (Mev) 


0. 0-1.0 
1 . 0 - 2. 0 

2. 0- 3.0 

3. 0- 4.0 

4. 0- 5.0 

5. 0- 6.0 

6. 0- 7.0 

7. 0- 8.0 


Np (Protons/cm 2 ) 


8.47E13 

1.23E11 

1.50E10 

7.99E9 

4.49E9 

2.50E9 

2.00E9 

1.10E9 

4.35E8 

4.49E8 
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8. 0- 9.0 

9.0- 10.0 

10.0-11.0 


3.98E8 











Table 2 (Continued) 


Energy Range Number 

AE (Mev) 

Np (Protons/cm 2 ) 

12 

11.0-12.0 

3.98E8 

13 

12.0-13.0 

3.98E8 

14 

13.0-14.0 

1.50E8 

15 

14.0-15.0 

1.14E8 

16 

15.0-30.0 

8.98E8 

17 

30.0-100.0 

8.47E8 

18 

>100.0 

1.40E8 


Table 3 

Solar-Flare Proton Flux Format 


Energy Range Number 

AE (Mev) 

Nfp (Protons/cm 2 ) 

1 

0. 0-1.0 

0.0 

2 

1. 0-2.0 

6.06E10 

3 

2. 0-3.0 

1.50E10 

4 

3. 0-4.0 

7.99E9 

5 

4. 0-5.0 

4.49E9 

6 

5. 0-6.0 

2.48E9 

7 

6. 0-7.0 

1.97E9 

8 

7. 0-8.0 

1.09E9 

9 

8. 0-9.0 

6.97E8 

10 

9.0-10.0 

4.49E8 

11 

10.0-11.0 

3.98E8 

12 

11.0-12.0 

3.98E8 

13 

12.0-13.0 

3.98E8 

14 

13.0-14.0 

2.99E8 

15 

14.0-15.0 

2.48E8 















Table 3 (Continued) 


Energy Range Number 

A E (Mev) 

Nfp (Protons/cm 2 ) 

16 

15.0-30.0 

1.79E9 

17 

30.0-100.0 

1.68E9 

18 

>100.0 

1.97E9 


Table 4 

Solar-Flare Alpha Particle Flux Format 


Energy Range Number 

A E (Mev) 

Nfa (Alpha-Particles/cm 2 ) 

1 

16-18 

3.98E7 

2 

18-20 

2.99E7 

3 

20-22 

2.99E7 

4 

22-25 

' 2.99E7 

5 

25-30 

3.98E7 

6 

30-32 

1.97E7 

7 

32-35 

1.97E7 

8 

35-40 

1.97E7 

9 

40-45 

1.97E7 

10 

45-47 

9.86E6 

11 

47-52 

1.17E7 

12 

52-57 

‘ 1.28E7 

13 

57-60 

. 4.96E6 

14 

60-80 

3.20E8 

15 

80-100 

1.97E7 

16 

100-200 

2.23E7 

17 

200-400 

4.49E5 

18 

>400 

9.86E5 














Table 5 

Damage Fa- tors for Electrons (KDE) 




























Table 6 (Continued) 


y Range 

Shielding Number 

mber 

0.5 

1.0 

1.5 

3.0 

6.0 

9.0 

200.0 

4 

5400.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

4500.0 

3000.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6 

3900.0 

3700.0 

2000.0 

0.0 

0.0 

0.0 

0.0 

7 

3500.0 

3500.0 

3100.0 

0.0 

0.0 

0.0 

0.0 

8 

3100.0 

3100.0 

3100.0 

200.0 

0.0 

0.0 

0.0 

9 

2800.0 

2800.0 

2800.0 

1400.0 

0.0 

0.0 

0.0 

10 

2700.0 

2700.0 

2700.0 

2000.0 

0.0 

0.0 

0.0 

11 

2600.0 

2600.0 

2600.0 

2100.0 

0.0 

0.0 

0.0 

12 

2500.0 

2500.0 

2500.0 

2100.0 

100.0 

0.0 

0.0 

13 

2500.0 

2500.0 

2500.0 

2100.0 

1000.0 

0.0 

0.0 

14 

2500.0 

2500.0 

2500.0 

2000.0 

1400.0 


0.0 

15 

2500.0 

2500.0 

2500.0 

2000.0 

1500.0 

100.0 

0.0 

16 

2500.0 

2500.0 

2500.0 

2000.0 

1800.0 

1500.0 

0.0 

17 

2300.0 

2300.0 

2300.0 

2000.0 

2000.0 

2000.0 

0.0 

18 

1000.0 

1000.0 

1000.0 

1000.0 

1000.0 

1000.0 

0.0 


This relationship states that an alpha particle of four times the energy of a pro- 
ton does 3.8 times the damage of the proton. Table 7 lists damage factors for 
alpha particles. The program user must specify an equivalent thickness of back 
shielding and coverglass for each computation; Table 8 lists the stored shielding 
densities, and equivalent thicknesses of fused silica coverglass and aluminum 
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Table 7 

Damage Factors for Alpha Particles (KDA) 


— 

Energy Range 
Number 

1 

Shielding Number 

0.0 

1.0 

1.5 

3.0 

6.0 

9.0 

200.0 

1 

1 ! 

20000.0 

10000.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

1 

15000.0 

13000.0 

0.0 

0.0 

0.0 

0.0 

0.0 

! 5 
3 ; 

i 

14500.0 1 

14000.0 

7000.0 

0.0 

0.0 

0.0 

0.0 

i 

4 ! 

13700.0 

13700.0 

11000.0 

0.0 

0.0 

0.0 

0.0 

5 

12500.0 

12500.0 

12000.0 

0.0 

0.0 

0.0 

0.0 

6 

11500.0 

11500.0 

11500.0 

2500.0 

0.0 

0.0 

0.0 

i 

7 

11000.0 

11000.0 

11000.0 

5200.0 

0.0 

0.0 

1 

o.o ; 

8 

10400.0 

10400.0 

10400.0 

7200.0 

0.0 

0.0 

! 

0.0 1 

9 

! ■ 

9800.0 

9800.0 

9800.0 

7900.0 

0.0 

0.0 

0.0 ' 

i 

1 10 ! 

9500.0 

I : 

9500.0 

9500.0 

00 

o 

o 

o 

o 

1700.0 

0.0 

0.0 | 

11 

| 

: 9500.0 

9500.0 

9500.0 

7800.0 

3650. 0 

0.0 

0.0 I 

12 

9500.0 

9500.0 

9500.0 

7700.0 

5300.0 

0.0 

0.0 

i I 

! 1 

! 13 

9500.0 

9500.0 

9500.0 

7600.0 

5700.0 

1600.0 

0.0 

1 14 

| 

9400.0 

9400.0 

9400.0 

7600.0 

6300.0 

4000.0 

0.0 

15 

9200.0 

9200.0 

9200.0 

7600.0 

7000.0 

6400.0 


16 

8600.0 

8600.0 

8600.0 

7700.0 

7400.0 

7100.0 

0.0 

17 

7000.0 

7000.0 

7000.0 

7000.0 

7000.0 

7000.0 

0.0 

18 

4000.0 

4000.0 

4000.0 

4000.0 

4000.0 

4000.0 

0.0 
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Table 8 

Shielding Numbers 


Shielding Density 
(gm/cm 2 ) 

0.016 

0.033 

0.05 

0.1 

0.2 

0.3 

* 

Equivalent mils 
of fused silica 
coverglass 

3.0 

6.0 

9.0 

18.0 

36.0 

54.0 

— 

Equivalent mils 
of Aluminum 

2.5 

6.0 

7.5 

15.0 

30.0 

45.0 

1000.0* 

Shielding Number 
for 

Computer Lookup 

0.5 

1.0 

1.5 

3.0 

6.0 

9.0 

200.0 


"'Simulates Infinite backshlelding 


backshielding. The user may specify any desired thickness of coverglasS or 
backshlelding within the range of stored data. The program will automatically 
make a linear interpolation between the damager factor values associated with the 
two closest thicknesses to determine the damage factor to be used. 

Once the charged particle population of each type of particle has been en- 
tered into the computer, and the coverglass and backshielding thicknesses have 
been selected, the computer can calculate the total equivalent 1-Mev electron 
flux using the relationship: 

<^ T =< 0 e + <Pp + <Pfp + <£f a (reference 1) 


where 

<P T = total 1-Mev equivalent flux 

<p e = electron 1-Mev equivalent flux 

Y>p = proton 1-Mev equivalent flux 

i/j fp = solar-flare proton 1-Mev equivalent flux 

<fi {a = solar-flare alpha particle 1-Mev equivalent flux 

The contribution of each particle to the total flux is calculated separately by a 
summation process. For example, <P p is calculated by multiplying the popu- 
lation of each of the 18 proton energy ranges by the appropriate damage factors 
for backshielding and coverglass, then summing each of these products. The 



following expression demonstrates this summation process: 


*p" £ [KDPXNP| cg + KDPXNP| bs ] 


Ae c i 


where 


KDP - proton damage factor fur a given energy range 

NP = number of protons in a particular energy range 

The equivalent flux for each of the other types of particles is calculated similarly. 
The value of <^> T is stored by the computer until called for by subroutine 
DEGRAD. 

C. Subroutine DEGRAD - Calculation of Solar Cell I-V Curve Irradiation 
Degradation 

Subroutine DEGRAD degrades a single solar cell I-V characteristic to 
compensate for damage caused by an equivalent 1-Mev flux. The input to this 
subroutine is either the output of PHI or a user specified 1-Mev flux; in the latter 
case subroutine PHI is not called and the computer immediately enters subroutine 
DEGRAD to degrade the I-V curve for the given flux. 

The undegraded solar cell 1-V curve is entered into the computer as a 
series of voltage and current points at room temperature and nominal solar in- 
tensity (140mw/cm 2 ). The computer interpolates the orginal points to arrive 
at 100 current and voltage pairs (Vi and Ii) at 10-mv increments. These inter- 
polated points will then, be shifted to account for the charged -particle irradiation 
degradation. 

The effect of a charged particle flux on the solar cell is a decrease in solar 
cell output current and voltage and a change in the I-V curve shape. Flight and 
laboratory experiments have shown that a solar cell will incur a ’’softening" of 
the I-V curve due to junction damage if low energy protons are allowed to reach 
the silicon where the coverglass does not completely protect the exposed cell 
surface. However, the computer program assumes an adequate shielding decign 
and therefore accounts for the effect of bulk damage only (decrease of minority 
carrier lifetime in the base material). The I-V curve shape change due to the 
bulk damage is a "squaring" of the curve shape. This change is implemented by 
a series resistance correction technique in the computer program. A very lucid 
graphical example of the curve shape change can be found in the Cherry and 
Statler report on U.S. and European solar cells (Reference 4). This reference 
includes I-V curves for Heliotek, Centralab and Texas Instruments 10-ohm-cm 
cells and Heliotek 2-ohm-cm cells, at 1-Mev electron flux levels from 0 to 10 16 
e/cm . When the 10 I-V curve is traced on a transparent paper overlay and 





this curve is now translated to the 0 flux curve, the sharper line of the irradiated 
I-V curve is quite obvious for till of the U.S. cells. 

Implementation of this effect in the computer is shown in Figure 2. 


I 



Figure 2. Curve Shape Correction 

A curve shape correction, dependent on the equivalent 1-Mev flux, effec- 
tively sharpens the knee of the curve by adding a negative "radiation resistance" 
(RR) correction to each voltage point of the curve: 

Vicorrected = Vij n jtj a | + Ii x RR (reference 1) 

RR is defined as "volts at short circuit current (VISC) divided by the initial, un- 
degraded value of short circuit current (ISCO)." Values of VISC, dependent 
upon flux level and solar cell base resistivity, are entered as a portion of the 
user specified input data (STINT table number 10: "NMVTBL CURVE SHAPE 
VOLTS VS PHI"). The proper value of VISC will be determined by a linear 
interpolation between the values associated with the flux levels closest to actual 
value of y9 T . Figure 2 shows the effect of the curve shape correction. Curve 
1 is the original I-V curve and curve 2 shows the sharpening of the knee caused 
by tlie charged-particle irradiation. 

Next the program shifts the total I-V curve parallel to the current axis an 
amount A I to account for the reduction in short circuit current caused by the 
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equivalent flux. The current increment is calculated as ISCO(l-Rl), where RI 
is defined as the relative short-circuit current degradation factor. Equations 
for RI, dependent on cell base resistivity and flux level, are contained in the 
program. Curve 3 of Figure 3 is determined by subtracting A I from each cur- 
rent point of curve 2: Ii — Ii-A I. The point PZCD is not the actual degraded 
open-circuit voltage, but merely a temporary value of the "point of zero current" 
of the cell. 


I 



A final characteristic correction is made to account for the reduction of 
open-circuit voltage in the solar cell caused by the equivalent flux. The relative 
open-circuit voltage degradation factor, RV, is calculated through the use of 
stored equations. (Unlike the previous degradation parameters, RV depends only 
on the equivalent flux and not on the solar cell base resistivity.) The degraded 
open-circuit Voltage (VOCR) is calculated as (PZCU x RV), where PZCU is the 
initial undegraded value of cell open-circuit voltage. The curve is then shifted 
parallel to the voltage axis an amount A V - defined as (PZCD-VOCR) - by sub- 
tracting A V from each voltage point: Vi - Vi- A V. Curve 4 of Figure 4 shows 
the fully degraded solar cell I-V curve. 
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Figure 5 shows typioal solar cell I-V curves for various flux levels . A 
small computer program employing the techniques of PHI and DEGRADE has been 
developed and is described in Reference 5. 

D . Subroutine STASH - Solar Cell Degradation and Temperature Effects 

The purpose of subroutine STASH is to account for all factors , with the 
exception of charged particle irradiation degradation, which effect the solar 
cell’s I-V characteristic. External series resistance effect, illumination in- 
tensity change, voltage measurement error and temperature effect are all in- 
cluded. The input to this subroutine, curve 1 of Figure 6, is either the output 
of DEGRAD or, if beginning-of-life calculations are being performed, the pro- 
gram’s original input solar cell I-V curve. In the latter case, subroutines PHI 
and DEGRAD are bypassed and the computer immediately enters subroutine 
STASH to account for the user-specified current and voltage degradation factors 
and temperature effects . All design uncertainties , or factors for a conservative 
design, should be included when the degradation factors and temperature coef- 
ficients are selected. 
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figure 4. Open Circuit Voltage Correction 



w 


A 





■ 



The I-V characteristic is first shifted parallel to the current axis an 
amount DELI to account for the current degradation factors Dll through DI4; each 
factor is defined as a percentage of short circuit current remaining after degra- 
dation. (For example, if it is assumed that Dll causes a 10% loss in short cir- 
cuit current, Dll would have a value of 90. 0.) DELI is defined as ISC 1(1- 7 ), 
where 7 is merely the product of the four current degradation factors . Although 
the four parameters usually refer to standard cell error, solar intensity variation, 
ultraviolet degradation and current measurement error, respectively, any four 
factors which effect the solar cell’s short circuit current output may be used. 

The current- shifted I-V curve is translated parallel to the voltage axis until V 2 
is intercepted, which accounts for the change in open-circuit voltage caused by 
the relative decrease in short circuit current. Curve 2 of Figure 6 shows the 
current -degraded I-V characteristic . 
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A voltage degradation factor which effects all voltages on the I-V curve 
equally is usually defined as a percent age of the open circuit voltage of the ref- 
erence I-V curve. Such a factor may be a voltage measurement error. The I-V 
curve is corrected to account for this factor by translating the entire curve par 

. t .1 _ » * 


allel to the voltage axis an amount DELV : DELV - 


percentage degradation 


100 
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Curve 3 of Figure 7 shows the voltage degraded I-V curve; VOC is the degraded 
value of open circuit voltage . 


A final correction is made to the cell’s I-V characteristic to account for 
an effective increase in the external series resistance. This increase, which 
is usually caused by either a series wiring loss or thermal cycling degradation, 
is defined as a percentage loss of voltage at the maximum power point of the 
solar cell I-V curve. There is, in addition, a corresponding voltage loss to all 
other points on the I-V curve which is proportional to the current at each point. 
Figure 8 shows the I-V characteristic before and after the resistance change; 
Vpmo JS th® original maximum power point voltage, while V pmr is the degraded 
value. 

After applying the degradation factors , subroutine STASH uses the current 
and voltage temperature coefficients to expand the fully degraded I-V curve into 
a family of curves at fifteen temperatures. Maximum power, maximum power 
point voltage a»d current, open circuit voltage and short circuit current values 
are calculated at each temperature . 

A detailed technical discussion of the techniques employed in this sub- 
routine is contained in Reference 6. 

E. Subroutine STINT - Data Storage and Retrieval 

Subroutine STINT merely stores in tables, and then supplies on demand, 
the values of variables which are functions of one, two or three arguments. The 
first use of STINT loads the tabular data into the computer. Subsequent STINT 
calls will ask for a linear interpolation to be performed which corresponds to 
the supplied argument values . 

As aa example, STINT table No. 12 of Appendix B gives a typical solar 
array temperature verses time profile for a Nimbus orbit. If the program user 
has requested that a solar array I-V curve be generated at 21 minutes into the 
sunlight portion of the orbit, MAIN will call STINT to determine the array tem- 
perature at this time in orbit. By interpolating between the array temperatures 
corresponding to the two closest time references - in this case 20.0 and 25.0 
minutes - STINT calculates an array temperature of 19.0 degrees; this value is 
returned to MAIN. 
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Figure 6. Short Circuit Current Degradation Effects 


I 



VOLTAGE 

Figure 7. Open Circuit Voltage Degradation Effects 
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Figure 8. External Series Resistance Correction 

One and two argument functions may be stored in single tables - Tables 
1 and 7 are examples of one and two argument functions, respectively. Three 
argument functions must be stored in consecutive two argument tables; the num- 
ber of two argument tables will be equal to the number of arguments values. 

Subroutine STINT is called by means of a statement with the following 
format: 

CALL STINT (ARG 1, ARG2, ARG3, FCT, KEY, NGRIPE, MINTBL, MAXTBL) 

The first three dummy variables are the three arguments; zeroes must be sub- 
situted for unused arguments. FCT is the variable value to be calculated, and 
KEY tells the computer what kind of operation is to be performed. A value of 
-1 indicates the table-loading mode, while a +1 calls for a linear interpolation 
to be performed. NGRIPE is an error flag. MINTBL and MAXTBL are the 
number of the STINT tables used; both tables are the same if the variable is a 
function of two arguments. MINTBL is set equal to the table number which con- 
tains the lowest value of arguments values and MAXTBL is set to the table which 
contains the highest argument 3 values if the variable is a function of three 
arguments . 
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HI. PROGRAM USAGE 

The following describes the mechanics of using the program in "non- 
programmer" language. Appendix A contains a program deck listing, while 
Appendix B illustrates a typical data deck set-up written in the FORTRAN IV 
computer language. 

The assembly of the complete program as it is submitted to the computer 
is shown in Figure 9. This assembly basically consists of two parts — a program 
deck and a data deck. The program deck is always used and is placed first in 
the assembly. It contains the MAIN routine and all of the subroutines used in 
the program (STINT, PHI, DEGRAD and STASH) and does not require any card 
changes to perform its function. 

The data deck contains all the numerical information the program requires 
for computation and defines the user-selected options for each rim. Consequently, 
the data deck must be prepared specifically for each run, or series of chained 
runs, to be made. Cards and tables in the data deck must be positioned in the 
order shown in the program assembly in Figure 9. The data deck description 
and format are presented below in the proper assembly sequence. 

Asterisk Data Card . The first card in the data deck is labeled: 

Column 1 * 


Columns 7-10 

DATA 

Date and List Option Card 

Column 1-2 

Number of month 

Column 3-4 

Number of day 

Column 5-6 

Number of year 

Column 9 

0 - Input data (STINT) tables not printed out 

Column 9 

1 - Input data (STINT) tables printed out in 
exponential format 
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blank carp 


NCOPES 
LAST RUN 


LABEL CARD 
FOR LAST RUN 
IN CHAIN 


BLANK CARD 


depending ON 

NUMBER OF PANELS 


SOLAR ARRAY 
DESCRIPTION 


END OF 
STINT TABLES 
CARD 


/ NCOOES 

FOR RUN NO. 1 

LABEL CARD 
FOR RUN NO. 1 


BLANK CARD 



^ STINT TABLE TITLE CARD I 


DATE CARD 
& NLIST OPTION 


PROGRAM DECK 


Figure 9. Program Assembly 


STINT Table Title Card 


Column 1 


Blank 


Column 2-72 Any Alpha-Numeric information 


Hi 
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The STINT tables are stacked one behind the other in the data deck in as- 
cending numerical order. The ten tables listed below must always be present 
in the data deck, in the order shown. 


Table # 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 


Table Name 
Solar Cell I-V Curve 

Relative Solar Array Current vs. Incidence Angle 

Orbital Electrons 1 yr. 600 N.M. (or suitable substitute) 

Orbital Protons 1 yr. 600 N.M. (or suitable substitute) 

Solar Flare Protons 1 yr. 600 N.M. (or suitable substitute) 

Solar Flare Alpha Particles 1 yr. 600 N.M. (or suitable 
substitute) 

Damage Factors for Electrons 
Damage Factors for Protons 
Damage Factors for Alpha Particles 
NMVTBL Curve Shape Volts vs . PHI 


The maximum number of STINT tables that the program will accept is 60; 
the 10 above tables plus, starting with Table No. 11, a sun angle vs. time table 
and a temperature vs . time table for each of the 25 panels comprising the solar 
array. If fewer than the maximum number (25) of panels is specified, fewer 
STINT tables are required. More than one panel may use the same angle or 
temperature table stored in STINT . 

The first card of each STINT table is a header card, which must be pre- 
pared in the following format: 


Columns 1-8 
Columns 9-12 

Columns 13-14 

Columns 15-16 

Columns 17-19 


Any alpha-numeric characters can be used for a date. 

Table number. Cannot be zero. Fixed point and 
right- justified. 

Number of argumenti values. Cannot be zero. Fixed 
point and right- justified. 

Number of argument 2 values. Cannot be zero, is 1 
for a function of one argument. Fixed point and 
right- justified . 

Not used . 
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Columns 20-70 : Any alpha-numeric characters desired. Usually used 

for table title. 

Columns 71-72 : 00 

After the header card, each card in the table uses 10 fields of 7 columns 
each for the argument values and the function values. The first card contains 
the first nine argument values in fields 2 through 10. In the following cards, 
field 1 contains an arguraent 2 value, and fields 2 through 10 contain correspond- 
ing function values. After all the argument 2 values have been spanned, the whole 
series of an argument) card followed by argument 2 cards can repeat until all 
the function values are used. If there is an argument^ value for the table, it 
goes into field 1 of the argument) card. Columns 71 and 72 on each card must 
contain a sequence number, starting with 01 for the first card. Figure 10 shows 
a typical STINT table coding sheet for a single argument (current as a function 
of voltage) STINT table. Figure 11 shows a two-argument STINT table format. 

After the last STINT table in the data deck, there is a card labeled END 
OF STINT TABLES, starting in Col. 21. Cols. 9-12 and 71-72 must be left 
blank on this card. This card is referred to as the blank card following the 
STINT tables in Figure 9. 

Run Label Card 

Following the END OF STINT TABLES card is a card containing any 
desired alpha-numeric information in columns 2-72, which usually describes the 
first rim to be made . 

NCODE 

Following the Run Label Card are the 24 NCODE cards. The card number, 
or NCODE, is right justified against Col. 3. The numerical value of the NCODE 
variable is left- justified against Col. 5 and must have a decimal point. Figure 
10 shows the NCODE names, the NCODE numbers, the NCODE values and a 
brief description of each NCODE for a sample computer run. Only the NCODE 
number and its numerical value are punched on the NCODE cards; the other 
data in Figure 12 are for information only. All 24 of the NCODES are initially 
loaded into memory, thus a single run or the first of a series of chained runs 
must contain all the NCODES in the data deck. The NCODE descriptions in 
Figure 12 are described below: 

a. TD is the total orbit suntime in minutes. 


b. DELTAT is the time increment in the orbit between solar array cal- 






















Figure 12. Format and Description of NCODES 



















calculating up to a maximum of 20 solar array I-V curves for each 
run. If the value of DELTAT is made larger than the value of TD 
(NCODE 1), only one solar array curve (at 0 minutes) will be computed. 

c. SIGISC is the solar cell short circuit current temperature coefficient. 

d . SIGVOC is the solar cell open-circuit voltage temperature coefficient 
and is punched on the card as a positive number? the program later 
gives it the proper negative sign. 

e. DI 1 is the first short circuit current degradation factor; it usually 
refers to a standard cell error. 

f. DI2 is the second short circuit current degradation factor; it usually 
refers to a solar illumination intensity variation. 

g. DI3 is the third short circuit current degradation factor; it usually 
refers to ultraviolet degradation. 

h. DI4 is the fourth short circuit current degradation factor; it usually 
refers to a current measurement error. 

i. DV 1 is the first maximum power point voltage degradation factor; it 
usually refers to an external series wiring loss . 

j. DV 2 is the second maximum power point degradation factor; it usually 
refers to a thermal cycling degradation loss. 

k. AVPMO is the maximum power point voltage of the input, undegraded 
solar cell. 

l. AIPMO is the maximum power point current of the input, undegraded 
solar cell. 

m. AVOCO is the open circuit voltage of the input, undegraded solar cell. 

n. THETA is the open circuit voltage degradation factor; it usually refers 
to a voltage measurement error . 

o. TNOT is the input solar cell reference temperature (Deg. C). 

p. ADIODE is the array blocking diode voltage drop, specified by user. 

A value of 0. 0 is punched if no blocking diode drop is desired in the 
output solar array curves . 
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q. DELTT is the temperature increment (Deg. C) between solar cell I-V 
curve calculations in STASH. 

r. ADDT is the temperature increment (Deg. C) to be added to TNOT to 
determine the highest STASH temperature . 

s. FLUX is punched as 0.0 if no charged particle irradiation damage is 
to be considered, in which case Subroutines PHI and DEGRAD do not 
compute. If FLUX has a value of 1-Mev equivalent flux (Example: 

2.6 x 10 14 is punched as 260000000000000.0 on the card) subroutine 
PHI is bypassed and subroutine DEGRAD will degrade the solar cell 
for the specified 1-Mev flux. If FLUX is given a value of -1.0, sub- 
routine PHI will compute a value of 1-Mev flux from tables of omni- 
directional particle fluxes that the user must supply in table locations 
3, 4, 5, and 6. This computed 1-Mev flux will then be automatically 
sent to subroutine DEGRAD for the appropriate sclar cell ,-rrections. 

t. CG is the sol ah cell coverglass thickness in mils of fused silica. 

u. BS is the equivalent solar cell backshielding thickness in mils of 
aluminum. 

v. BOHMS is the solar cell base resistivity (ohm-cm). 

w. NDGRAD must be set to 1. 0 in the first run. This causes the machine 
to automatically degrade the solar cell and expand it for temperature 
in subroutine STASH as specified by the degradation and temperature 
parameters in the NCODES. When chaining additional runs, if the 
solar cell degradations are not changed, NDGRAD should be set to 0.0 
in the second run . By setting NDGRAD to 0.0, needless repetitive 
computations in the solar cell subroutines are eliminated. 

x. NEND must be set to either 1.0, 2.0, 3.0, or 4.0 for each run, and 
determines which of the following options is selected: 

NEND = 1.0 Do solar cell calculations only (PHI, DEGRAD and STASH) 
and stop. 

NEND = 2.0 Do solar cell calculations and read new set of run 
instructions. 

NEND = 3.0 Do solar cell and solar array calculations and stop. 

NEND = 4.0 Do solar cell and solar array calculations and read new 
set of run instructions . 
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If no solar array calculations are to be made (NEND = 1. 0 or 2. 0), a blank 
card must follow NCODE 24. This tells the computer to stop reading in data and 
to start computing. If no additional (chained) run is to be made, this blank card 
is the last card in the data deck. 

Array Signal Card 


If NCODE 24 has been punched with a 3. 0 or a 4. 0, immediately following 
the NCODE 24 card must be a card containing 999 in columns 1-3. This card 
tells the computer that solar array information is to follow. 

NPANEL Card 


Following the Array Signal Card is the NPANEL card, which contains the 
number (NPANEL) of solar panels in the array (maximum number of panels is 
25). This number must appear right- justified in columns 1-3; no decimal point 
is required. 

Panel Description Cards 


Following the NPANEL card is a panel description card for each solar 
panel in the array, up to a maximum of 25 panels. The number of these panel 
description cards must agree with the value of NPANEL. Each card contains 
four fields of ten columns each, in floating point format (requires decimal point). 


Columns 

Variable 

Typical Value 

1-10 

No. of Series Solar Cells per String 

100.0 

11-20 

No. of Parallel Strings per Panel 

10.0 

21-30 

Panel Incidence Angle vs. Time 
Table Location in STINT 

11.0 

31-40 

Panel Temperature vs. Time 
Table Location in STINT 

12.0 


Following the Panel Description Cards is a blank card. This tells the 
computer to stop reading in data and to start computing. If it is desired to chain 
an additional run, a new Run Label Card and only those NCODES and Panel 
Description Cards that contain changed or new information should be placed after 
the blank card. As many runs as arc desired can be chained in this manner, en- 
suring that each new run starts with a Run Label Card and ends with a blank card. 
Refer again to Figure 9 for the proper sequence of card positions for chained runs . 
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IV. PROGRAM OUTPUT 


The information that the computer prints out after a run consists of the 
following items: 

1 . Stint Table Listing Option : There is an option in the program to list 
all the STINT tables. This is done by taking the date card at the be- 
ginning of the data package and either punching a one (1) or a zero (0) 
in column 9. Punching a one (1) will produce a listing of the STINT 
tables . A zero (0) will not produce a listing of the STINT tables . 

2. Input Data Page : 

a. Run number and date 

b. Run comments (as specified on input cards) 

c . Listing of NCODE numbers , names and values 

d. Solar array description: panel number, number of series solar 
cells, number of parallel solar cell strings, number of STINT table 
which contains the solar illumination incidence angle, number of 
STINT table which contains the array temperature-time profile 

3 . Subroutine PHI : 

a. The computed electron flux 

b. The computed proton flux 

c . The computed solar flare proton flux 

d . The computed solar flare alpha particle flux 

e. The computed total flux 

The fluxes are in equivalent 1-Mev electrons/cm 2 

4. Subroutine DEGRAD 

a. The solar cell I-V curve (irradiation degraded) 

b. Short circuit current 




s ; 


m 
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c . Current at the maximum power point 

d. Voltage at the maximum power point 

e. Open circuit voltage 

5. Subroutine STASH : Values of temperature for which the degraded solar 
cell I-V curve has been prepared are listed in a row across the page . 
Values of maximum-power voltage and current, open circuit voltage 
and short-circuit current appear in columns under each temperature . 
Values of every other calculated current and voltage pair comprising 
tiie I-V curve and stored in the computer memory, are listed in col- 
umns under each temperature (in 20 mv increments) • 

6. Solar Array Data : 

a. Value of blocking diode voltage drop. 

b. Time in orbit sunlight (up to 20 time increments starting with 0.0 
min.). 

c. Solar array voltage in one-volt steps from 0 volts to open-circuit 
voltage. 

d . Total solar array current at each voltage above for each time in- 
crement; a maximum of 20 solar array I-V curves can be produced 
per computer run. 


V. SUMMARY 

A computer program has been described which accounts for all factors that 
ali'eet the power output of a solar cell. Charged particle irradiation degradation, 
external series resistance, illumination intensity variation and temperature ef- 
fects are taken into account. Using the fully degraded solar cell characteristic 
as a basis, the program will determine the total I-V curve of solar arrays 
consisting of up to 25 separate panels. Each panel may have its own series- 
parallel solar cell arrangement, temperature versus time profile and solar 
incidence angle versus time profile. 
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GLOSSARY 

• DELI - decrease in short circuit current caused by current degradation 

factors; used in STASH. 

• DELV - decrease in open circuit voltage caused by the voltage degradation 

factor; used in STASH. 

• Ii - any particular current value of input I-V curve . 

• ISCO - initial, undegraded value of short curcuit current. 

• I, - radiation degraded value of short circuit current; Used in STASH. 

• I 2 - degraded value of short circuit current; used in STASH. 

• KD - arbitrary damage factor 
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• KDA - alpha particle damage factor; used in PHI. 

• KDE - electron damage factor; used in PHI. 

• KDP - proton damage factor; used in PHI. 

• Ne - number of electrons in a given energy range; used in PHI. 

• Nfa - number of solar-flare alpha particles in a given energy range; used 

in PHI. 

• Nfp - number of solar-flare protons in a given energy range; used in PHI. 

• Np - number of protons in a given energy range; used in PHI. 

• PZCD - degraded point of zero current; a temporary parameter used in 

DEGRAD. 

• PZCU - undegraded point of zero current; the initial value of open-circuit 

voltage used in DEGRAD . 

• RI - relative short circuit current degradation factor; used in DEGRAD. 

• RR - radiation resistance, defined as VISC/ISCO; used in DEGRAD. 

• RV - relative open circuit voltage degradation factor; used in DEGRAD. 

• Vi - any particular voltage point on input I-V curve. 

• VISC - volts at short circuit current; used in DEGRAD. 

• VOCR - degraded value of open circuit voltage; used in DEGRAD. 

• VOC - degraded value of open circuit voltage; used in STASH. 

• Vp mo “ initial value of maximum power point voltage; used in STASH. 

• Vpmr - the degraded value of maximum power point voltage; used in STASH. 

• Vj - radiation degraded value of open circuit voltage; used in STASH. 

• V, - degraded value of open circuit voltage; used in STASH. 


33 


• A I - decrease in short circuit current caused by the equivalent flux, de- 

fined as ISCO-ISCR; used in DEGRAD. 

• <p e - electron 1-Mev equivalent flux. 

• - solar-flare alpha particle 1-Mev equivalent flux. 

• ¥>fp - solar-flare proton 1-Mev equivalent flux. 

4 </>p - proton 1-Mev equivalent flux. 

• 'Pj - total 1-Mev equivalent flux. 
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»*** SOLAR ARRAY SYNTHESIS PROGRAM DEVELOPED 
***. PROGRAMMERS R. RASMUSSEN, P.J. HYLAND, 

DIMENSION NDATE (3) 

DIMENSION TC 20 ) 

DIMENSION XI A(20,100) 

DIMENSION XC(25),XS(25),NANGLE<25),NTEMPS<25> 


***+ READ DATE 
I A=5 
I B = 6 

READ ( IA,7000)NDATE( 
CALL SLITE (0) 

NRUN-1 

DO 6001 1*1,25 
XC(l)*0,0 
XS( I )»0.0 
NANGLE ( I ) *0 
6001 NTEMPSd )-0 
TD*0.0 

DELTAT-1000.0 

SIGISC*0.0 

SIGVOC-O.O 

011 * 100.0 

Dl 2*100.0 

Dl 3*100.0 

014*100.0 


, NDATE ( 2 ) , NDATE ( 3 ) , NL I ST 


DV1*100 .0 
DV2*100 .0 
AVPMO-O.O 
Al PMO“0 . 0 
AVOCO*0.0 
THETA*100.0 
TNOT-O.O 
ADIODE«0.0 
DELTT*0.0 
ADDT=0.0 
FLUX=0.0 
CG*0.0 
BS*0.0 
BOHMS*0.0 
NDGRAD-0 
NEND-0 
NCELLT-1 
NKOST-2 
NELECT-3 
NPROTT-4 
NFLARE-S 
NALPHA-6 
NKDE“7 
NKDP-8 
NKDA*9 
NMVTBL-10 
NKDFL«8 
KEY*NLI ST-1 

C* <*** READ TABLES COMMENT CARD 

100 READ ( I A, 7001) 

C**»* LOAD STINT TABLES 

CALL STINT (0.0, 0.0, 0.0, 0.0, 

IF (NGRIPE) 101,101,102 

102 WRITE (IB, 7001) 

WRITE (IB, 7002) 

PRINT 7002 
CALL EXIT 

(>*•* READ RUN COMMENTS CARD 

101 READ ( I A, 7003) 

C***» PRINT HEADER CARD 

WRITE ( IB, 700DNRUN, (NDATEC J), J - l»3) 


FOR NASA (GSFC) BY RCA-AED01 00010 
JUNE 1968 01 00020 

01 00030 
01 00040 
01 00050 
01 00060 
01 00070 
01 00080 
01 00090 
01 00100 
01 00110 
01 00120 
01 00130 
01 00140 
01 00150 
01 00160 
01 00170 
01 00180 
01 00190 
01 00200 
01 00210 
01 00220 
01 002 JO 
01 00240 
01 00250 
01 00260 
01 00270 
01 00280 
01 00290 
01 00300 
01 00310 
01 00320 
01 00330 
01 00340 
01 00350 
01 00360 
01 00370 
01 00380 
01 00390 
01 00400 
01 00410 
01 00420 
01 00430 
01 00440 
01 00450 
01 00460 
01 00470 
01 00480 
01 00490 
01 00500 
01 00510 
01 00520 
01 00530 
01 00540 
01 00550 
01 00560 
01 00570 
01 00580 
01 00590 
01 00600 
01 00610 
01 00620 
01 00630 
01 00640 
01 00650 
01 00660 


KEY, NGRIPE, 0,0) 


A-l 









WRITE (IB, 7001) 


01 


WRITE (IB, 7005) 


01 


WRITE (IB, 7003) 


01 


WRITE (IB, 7006) 


01 

c**** 

READ INPUT DATA,NCQDES 


01 

103 

READ ( IA,7O07)NCODE,PARAM 


01 


IF (NCODE-999) 1103,106,1103 


01 

1103 

IF (NCODE) 105,105,5000 


01 

106 

READ ( IA,7007)NPANEL 


o:i 


DO 6000 K«1,NPANEL 


01 


READ ( IA,6600)P1,P2,P3,P4 


01 

6600 

FORMAT ( 1QF10 . 5 ) 


01 


X C ( K ) n P 1 


01 


XS ( K) a P2 


ox 


NANGLE F K) “P3* .01 


01 


NTEMPS ( K ) n P4+ . 01 


01 

6000 

WRITE ( IB,6601)K,XC(K),XS(K),NANGLE(K),NTEMPS(K) 

01 

6601 

FORMAT ( 7H0PANEL 13/ 16H SERIES 

CELLS - F5 . 1, 20H PARALLEL 

STRINGS01 

1 « F6 ♦ 2/ 15H ANGLE TABLE * I2,14H 

TEMP TABLE - 12) 

01 


WRITE ( I B, 6602 ) 


01 

6602 

FORMAT ( 1H0) 


01 


GO TO 103 


01 

5000 

GO TO (1,2,3,4,5,6,7,8,9,10.11.12 

,13,14,15,16,17,18,19,20, 

01 

121, 22,23,24,25), NCODE 


01 

1 

TD-PARAM 


01 


WRITE ( IB,201)NCODE,TO 


01 


GO TO 103 


01 

2 

DELTAT-PARAM 


01 


WRITE ( I B, 202 ) NCODE, DELTAT 


01 


GO TO 103 


01 

3 

SIGISOPARAM 


4 01 


WRITE ( IB,203)NCODE,SIGISC 


01 


GO TO 103 


01 

4 

SIGVOOPARAM 


01 


WRITE ( 1 B, 204 )NCODE, $1 GVOC 


01 


GO TO 103 


01 

5 

DI1«PARAM 


01 


WRITE ( IB,205)NCODE,OI 1 


01 


GO TO 103 


01 

6 

D I 2«PARAM 


01 


WRITE ( IB,206)NCODE,OI2 


01 


GO \C 103 


01 

7 

DI3«PARAM 


01 


WRITE ( 1 B, 207 )NCODE, 01 3 


01 


GO TO 103 


01 

8 

DI4«PARAM 


01 


WRITE ( IB, 208) NCODE, DI4 


01 


GO TO 103 


01 

9 

DV1«PARAM 


01 


WRITE ( 1 B, 209 )NCODE, DV1 


01 


GO TO 103 


01 

10 

DV2«PARAM 


01 


WRITE ( IB,210)NCODE,DV2 


01 


GO TO 103 


01 

11 

AVPMO*PARAM 


01 


WRITE ( IB,211)NCOOE,AVPMO 


01 


GO TO 103 


01 

12 

Al PMO-PARAM 


01 


WRITE ( IB,212)NCODE,AIPMO 


01 


GO TO 103 


01 

13 

AVOCO-PARAM 


01 


WRITE ( !B,213)NCGDE,AVOCO 


01 


VOCO 1 «AVOCO 


01 


GO TO 103 


01 

14 

THETA-PARAM 


01 


WRITE ( IB, 214)NCOOE, THETA 


01 


00 6 7 0 

oor.80 
on ono 
00700 
00 7 10 
00720 
QQ730 
00740 
00700 
00700 
00770 
00780 
00 7 00 
00800 
00 810 
00870 
00830 
00840 

o n h & o 

008 GO 
00870 
00880 
00800 
00900 
00910 
00920 
00930 
00940 
00950 
00960 
00970 
00980 
00990 
01000 
01010 
01020 
01030 
01040 
01050 
01060 
C1070 
01080 
01090 
01100 
OHIO 
01120 
01130 
01140 
01150 
01160 
01170 
03180 
01190 
01200 
01210 
01220 
01230 
01240 
01250 
012G0 
01270 
01280 
01290 
01300 
01310 
01320 


THETA-THETA/100,0 
GO TO 103 

15 TNOT-PARAM 

WRITE <IB,215)NCODE,TNCT 
GO TO 103 

16 AD I ODE-PARAM 

WRITE (IB,216)NCODE,AD!OPE 
GO TO 103 

17 OEETT « PA RAM 

WRITE (IB,217)NCODE,DELTT 
GO TO 103 

18 ADDT-PARAM 

WRITE ( IB,218)NCODE,AODT 
GO TO 103 

19 FLUX-PARAM 

WRITE ( I B, 219 )NCODE, FLUX 
GO TO 103 

20 CG-PARAM 

WRITE ( I B, 220 )NCODE, CG 
GO TO 103 

21 BS-PARAM 

WRITE (IB,221)NCODE,BS 
GO TO 103 

22 BOHMS-PARAM 

WRITE ( I B, 222 )NCODE, BOHMS 
GO TO 103 

23 NOG RAD® PA RANH *01 

WRITE ( IB,223)NCODE,NDGRAD 
GO TO 103 

24 NEND-PARAM+,01 

WRITE (1B,224)NC0DE,NEND 
GO TO 103 

25 CONTINUE 
GO TO 103 

C **** END DATA LOADER 
O*** HOUSEKEEPING FOLLOWS 
105 IF (NDGRAO) 70,501,70 

70 Dt SC**(D 1 1/ 100 , )* (DI 2/100 , )*(DI 3/100 ♦ )*(D I 4/100 ♦ ) 

DV «OVl/ 100 . )* (DV2/ 100 • ) 

I F ( BOHMS- 5 • 0)81,81,82 

81 DENFAC a 0 • 065 
GO TO 83 

82 DENFAC-0* 0 

83 PHITOT-FLUX 

I F (FLUX ) 71, 73, 72 

71 CALL PHI (CG, BS,NKDE,NKDP,NKDFL,NELECT,NPROTT>NFLARE, 

INALPHA, NKDA, PHI TOT) 

FLUX-PH I TOT 

72 CALL DEGRAD ( PH I TOT, BOHMS, NCELLT,HMVTBL) 

CALL SLITE (1) 

C**** INITIAL ENTRY INTO SUBROUTINE STASH, I NITI ALI ZATI ON 

73 CALL STASH (D I SC, THETA, AVPMO, Al PMO, TNOT, -1 ,NGR I PE,NCELLT, VOCO I , 
1ADDT,DELTT> 


01 01330 
01 01340 
01 01350 
01 01360 
01 01370 
01 01380 
01 01390 
01 01400 
01 01410 
01 01420 
01 01430 
01 01440 
01 01450 
01 01460 
01 01470 
01 01480 
01 01490 
01 01500 
01 01510 
01 01520 
01 01530 
01 01540 
01 01550 
01 01560 
01 01570 
01 01580 
01 01590 
01 01600 
01 01610 
01 01620 
01 01630 
01 01640 
01 01650 
01 01660 
01 01670 
01 01680 
01 01690 
01 01700 
Cl 01710 
01 01720 
01 01730 
01 01740 
01 01750 
01 01760 
01 01770 
01 01780 
01 01790 
01 01800 
01 01810 
01 01820 
01 01830 
01 01840 
01 01850 



SECOND ENTRY , M f 0 SUBROUTINE STASH, 1 N 1 T 1 AL I ZATI ON 

01 

01860 


CALL STASH (DV, Dl SC, S 1 G 1 SC, S 1 GVOC, TNOT, 0, NGR 1 PE, NCELLT, 

DENFAC, ADDT01 

01870 


1, DELTT) 

01 

01880 


1 F (NGRI PE)74, 76,74 

01 

01890 

74 

WRITE (IB, 7009) 

01 

01900 


CALL EXIT 

01 

01910 

76 

GO TO (80, 30, 501, 501 ), N END 

01 

01920 

30 

NRUN*NRUN+1 

01 

01930 


GO TO 1C1 

01 

01940 

80 

CONTINUE 

01 

01950 


CALL EXIT 

01 

01960 


A-3 


201 FORMAT (17,50m TO, ORBIT SUNTIME DURATION (MINUTES) 01 

1F12.4) 01 

202 FORMAT (I7,5GH DELTAT, T (ME BETWEEN CALCULATIONS (MINUTES) 01 

IF 12 * 4 ) 01 

203 FORMAT (I7,5GH $IGISC,ISC TEMP . COEF . (AMPS/DEG , C) 01 

IFI2.6) 01 

204 FORMAT (17,50m SIGVOC,VOC TEMP* COEF • (VOLTS/ DEG * C ) 01 

1F12.6) 01 

205 FORMAT (I7,5QH Dll, ISC DEGRADATION FACTOR (PERCENT) 01 

1 F 1 2 • 4 ) 01 

206 FORMAT (I7,50H DI2,ISC DEGRADATION FACTOR (PERCENT) 01 

1 F12 ♦ 4 ) 01 

207 FORMAT (I7,50H DI3,!SC DEGRADATION FACTOR (PERCENT) 01 

1F12 * 4 ) 01 

208 FORMAT (!7,50M DI4,ISC DEGRADATION FACTOR (PERCENT) 01 

IF 12 * 4 ) 01 

209 FORMAT C I 7, 50H DV1,VPM DEGRADATION FACTOR (PERCENT) 01 

1 F 1 2 4 4 ) 01 

210 FORMAT (I7,50H DV2,VPM DEGRADATION FACTOR (PERCENT) 01 

1 F 1 2 4 4 ) 01 

211 FORMAT (17, 50H AVPMO,VPM OF UNDEGRADED CELL (VOLTS) 01 

1F12 . 4 ) 01 

212 FORMAT (I7,50H AIPMO,IPM OF UNDEGRADED CELL (VOLTS) 01 

l c 12 , 4 ) 01 

213 FORMAT (I7,50H AVOCO,VOC OF UNDEGRADED CELL (VOLTS) 01 

1F12 * 4 ) 01 

214 FORMAT (I7,5UH THETA, VOC DEGRADATION FACTOR (PERCENT) 01 

1F12 4 4 ) ' 01 

215 FORMAT (I7,50H TNOT, SOLAR CELL REFERENCE TEMP* (DEG . C) 01 

1F12 * 4 ) 01 

216 FORMAT (!7,50H ADIODE, ARRAY BLOCKING DIODE DROP (VOLTS) 01 

1F12 * 4 ) 01 

217 FORMAT (I7,50H OELTT , TEMP. INCREMENT FOR STASH STORAGE (DEG* C) 01 

1F12.4) 01 

218 FORMAT (!7,50H ADDT, HIGHEST STASH TEMP* MINUS TNOT (DEG.C) 01 

1F12 . 4 ) 01 

219 FORMAT (I7,50H FLUX, DENI 1-MEV ELECTRON FLUX ( ELEC/SQ* CM* ) 01 

1E12 * 4 ) 01 

220 FORMAT (I7,50H CG, COVER GLASS THICKNESS (MILS) 01 

1F12.4) 01 

221 FORMAT (I7,50H BS, BACKSH I ELD I NG THICKNESS (MILS. OF ALUM.) 01 

If 12 . 4 ) 01 

222 FORMAT (I7,50H BOHMS, SOLAR CELL BASE RESISTIVITY (OHM-CM) 01 

1F12.4) 01 

223 FORMAT (I7,50H NDGRAD, INITIALIZE STASH (1.0 OR 0.0) 01 

118) 01 

224 FORMAT (I7,50H NENO, END OF RUNS KEY (1.0, 2. 0,3.0 OR 4.0) 01 

118) 01 

225 FORMAT (I7,50H NLIST, STINT TABLE LISTING OPTION (1.0 OR 0.0) 01 

118) 01 

C**** FORMATS 01 

7000 FORMAT (3I2,2X,I1) 01 

7001 FORMAT ( 72H 01 

1 ) 01 

7002 FORMAT (34H UNABLE TO READ TABLES, ABORT RUN) 01 

7003 FORMAT (72H 01 

1 ) 01 

7004 FORMAT ( 12H1 RUN NO.I3,18H ON THIS DATE OF 12, 1H- 12, 1H- 1 2, 75H 01 

1 SOLAR ARRAY SYNTHESIS PROGRAM 01 

2 / 20H0 TABLE COMMENTS.) 01 

7005 FORMAT (18H0 RUN COMMENTS.) 01 

7006 FORMAT (30H0 NEW OR CHANGED PARAMETERS/1H0) 01 

7007 FORMAT ( I 3, IX, F20 . 0 ) 01 

7008 FORMAT (35H0 ERROR IN TABLE LOOKUP, NGRIPE «tl) 01 

7009 FORMAT (25H STASH DID NOT INITIALIZE) 01 

O*** START SOLAR ARRAY COMPUTATIONS 01 


01070 
01080 
01900 
02000 
02010 
02020 
020 50 
02040 
02050 
02060 
02070 
02080 
02090 
02100 
02110 
02120 
02130 
02140 
02150 
02160 
02170 
02180 
02190 
02200 
02210 
02220 
02230 
02240 
02250 
02260 
02270 
02280 
02290 
02300 
02310 
02320 
02330 
02340 
02350 
02360 
02370 
02380 
02300 
02400 
02410 
02420 
02430 
02440 
02450 
02460 
02470 
02480 
02490 
025C0 
02510 
02520 
02530 
02540 
02550 
02560 
02570 
02580 
02590 
02600 
02610 
02620 


A-4 


501 DO 503 INCR»1,20 
DO 503 JVOLT-1,100 

X I A( INCR,JVOLT) - 0.0 
503 CONTINUE 
T( 1) - 0.0 
(NCR - 1 
GO TO 599 

502 T ( I NCR+1 ) « T ( I NCR) ♦ DELTAT 

IF (TD - T(INCR*1)) 700,505,505 
505 (NCR * I NCR ♦ 1 

IF (20 - I NCR) 700,599,599 
C K IS PANEL NUMBER 

599 K*1 

600 CALL STINT (T ( I NCR) , 0 . 0, 0 .0, PANANG, 1, NGR I PE,NANGLE< K) , NANGLE ( K) ) 
IF (NGRI PE )690, 601,690 

690 NGRIPE-1 
GO TO 1000 

601 CALL STINT (T( (NCR) ,0 .0,0 .0, PTEMP, 1,NGRI PE,NT£MPS(K) ,NTEMPS( K) ) 
IF (NGRI PE) 691, 602, 691 

691 NGRIPE-2 
GO TO 1000 

602 CALL ST ! NT ( PANANG, 0 .0,0 ♦ 0, EFFECT,!, NGRI PE, 2, 2 ) 

IF (NGR I PE ) 692, 603, 692 

692 NGR I PE-3 
GO TO 1000 

603 JVOLT-1 

IF (EFrECT) 607,607,604 

604 AVOLT-FLOAT (UVOLT)-l . 0+AD IODE 
ARG1-AVOLT/XC (K) 

I F (JVOLT-1 ) 610, 610, 611 

C**** GET VALUE OF FULL INTENSITY ISC FROM STASH 

610 CALL STASH (TLU, 0.0, 0.0/ PTEMP,0 .0, 1, 0 .0,0.0, 0*0, 0.0, DELTT) 
DELJ«TLU*(1. 0-EFFECT) 

PHIL-1. 0/EFFECT 

DELVOC»(0.026*(273.0+PTEMP)/300.0)*ALOG (PHIL) 

C**** GET VALUE OF FULL INTENSITY VOC FROM STASH 

CALL STASH(OCV,0.0,0.0,PTEMP, 0.0, 989, 0.0, 0.0, 0.0, 0.0, DELTT) 

C**** GET VALUE OF VOLTAGE AT CURRENT DELI FROM STASH 

CALL STASH(VAI,0.0,DELI,PTEMP, 0.0, 988, 0.0, 0.0, 0.0, 0.0, DELTT) 
DELTAV-OCV-VA I -DELVOC 

611 I F( ARG1-DELTAV) 606, 605/ 605 

605 ARG1-ARG1-0ELTAV 

606 CALL STASH(CELL I, 0 , 0, ARG1, PTEMP, 0.0, 1,0, 0,0. 0,0. 0,0.0, DELTT) 
PANEL I -X$(K)* (CELL! -DEL! ) 

I F (PANEL I ) 607, 607, 609 

607 PANEL 1-0.0 

608 K-K+l 

I F (NPANEL-K) 502, 600, 600 

609 XIA(INCR,JVOLT)-XIA(INCR,JVOLT)*PANELI 
JVOLT-JVOLT+1 
!F(100-JVOLT)608,604,604 

1000 WRITE ( I B, 7008 )NGRI PE 
CALL EXIT 

700 WRITE (IB,9000)ADIODE 
T ( I NCR+1 )-0. 0 
I BEG-1 
I F.ND-10 
DO 400 K-l, 2 
AVOLT-O. 0 

WRITE ( 18,9001) (T( I NCR) , I NCR-1 BEG, I END) 

401 DO 410 J-1,100 

WRITE ( I B, 9002 ) AVOLT, (X I A( I NCR, J ) , I NCR- 1 BEG, I END) 

AV0LT-AV0LT*1.0 

402 DO 420 L-l BEG, I END 

IF (XIA(L,J)-0. 0)410, 420,410 
420 CONTINUE 


01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 


02740 

02750 

02760 

02770 

02780 

02790 

02800 

02810 

02820 

02830 

02840 

02850 

02860 

02870 

02880 

02890 

02900 

02910 

02920 

02930 

02940 

02950 

02960 

02970 

02980 

02990 

03000 

03010 

03020 

03030 

03040 

03050 

03060 

03070 

03080 

03090 

03100 

03110 

03120 

03130 

03140 

03150 

03160 

03170 

03180 

03190 

03200 

03210 

03220 

03230 

03240 

03250 

03260 

03270 

03280 


GO TO 430 
410 CONTINUE 
430 IBEG-11 
I END-20 

WRITE (16,9005) 

400 CONTINUE 
NRUN-NRUN+ 1 

GO TO (80,80,80,101),NEND 

9000 FORMAT (60H1 SOLAR ARRAY l-V CURVES WI01 

1TH <F4.2,28H) VOLTS BLOCKING DIOOE DROP ) 01 

9001_FORMAT_(/17H TIME IN SUNLIGHT/ 11H (Ml MUTES) 10F10 . 2, //51X, 16HCURRE01 


01 

01 

01 

01 

01 

01 

01 

01 


1NT( AMPERES )/2X,5HVOLTS) 

9002 FORMAT ( FS . 1, 6X, 10F10 . 2 ) 

9003 FORMAT (1H1) 

END 

SUBROUTINE STINT ( ARG1, ARG2, ARG3, FCT, KEY,NGR I PE, Ml HTBL,MAXTBL) 
DIMENSION NUMPT$( 61 ),L1(60),L2(60),L3(60), STG( 3600), DUMMY (10) 
DIMENSION NAME (14) 

EQUIVALENCE (NAT,L3(1)) 

C**«* SIZE OF STG IS CALCULATED BY SUM OF ((1+N(ARG1))*(1*N(ARG2))) 

I A* 5 
IB*6 

NGR I PE*0 
IF (KEY) 1,1,70 
646 FORMAT (1X,10E11 # 4, !2) 

1 NG*1 
NORMAL*! 

WRITE (18,1357) 

1357 FORMAT (34H0 TABLE DATE CONTENTS) 

GO TO 55 
2000 NG*2 

NORMAL *2 
WRITE (IB, 1257) 

1257 FORMAT (1H1) 

3000 RETURN 

775 NGR I PE*1 
RETURN 

776 NGR I PE*2 

WRITE ( I B, 9000 )ARG1, ARG2, ARG3,MI NTBL,MAXTBL 
RETURN 


18H MINTBL*I4,8H MAXTBL* 14) 

; BEGINNING OF STINT 
55 NUMTBL*1 
NUMPTS ( 1)^0 

102 READ (IA,57)DA1, DA2, DA3, K, LI (NUMTBL) , L2 (NUMTBL) ,NAME, 

57 FORMAT (A2, A3, A3, I 4, 2 I 2, 13A4, A2, 1 2 ) 

WRITE ( IB, 1157) K, DAI, DA2,DA3, NAME 

1157 FORMAT (/ I 8,5X,A2, A3, A3, 5X, 13A4, A2) 

104 IF(ISEQ) 69,58,69 

58 IF (K) 99,99,1159 

1159 IF (K-60) 59,59,1103 

59 L8*L1 (NUMTBL ) 

Nl*(L8-l)/94l 
DO 68 I S*l, N1 
NAT* ( I S-l )*9+l 

IF (IS-N1) 60,61,60 

60 L4*NAT*8 
GO 70 62 

61 14*18 

62 L5»NUMPTS(NUMTBL)*1 
L6*15*NAT 
L7*L5+L4 

JJ«0 

L9*L2 ( NUMTBL ) 

LM*L5*L8 


I SEQ 


01 

01 

01 

01 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 


9000 FORMAT (20H0 ERROR IN TLU, ARG1*E12. 5, 6H ARG2*E12*5,6H ARG3*E12.5, 02 

lOIJ imiTDI.il. an w 


02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 


03290 

03300 

03310 

03320 

03330 

03340 

03350 

03360 

03370 

03380 

03390 

03400 

03410 

03420 

03430 

00020 

00030 

00040 

00050 

00060 

00070 

00080 

00090 

00100 

00110 

00120 

00130 

00140 

00150 

00160 

00170 

00180 

00190 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00560 

00370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 

00490 

00500 

00510 

00520 
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LN-UHL9 

105 READ ( IA,64)(DUMMY(K),K«1,10), ISEQ 

64 FORMAT (10E7.0, 12) 

IF (KEY) 107,11,11 

11 WRITE (IB,646)( DUMMY (K),K-1,10), ISEQ 

107 STG ( L5 ) a DUMMY ( 1 ) 

K-2 

DO 65 J»l6,t7 
STG ( J ) *DUMMY ( K) 

65 K-K+ 1 

IF (ISEQ-(US-1)*(L9*1)4JJ*1)) 69,66,69 

66 t6 c LN*NAT 
L7-LN+ 14 
L5-IM+1+JJ 

IF (JJ-L9) 67,68,69 

67 JJ-JJ+1 
LN-LN+L8 
GO TO 105 

68 CONTINUE 

109 LEE-NUMPTS (NUMTBl )♦ ( L8+1)*(L9+1) 

IF UEE-3600) 1100,1100,1101 

1100 IF (NUMTBL-60) 1102,108,1103 

1102 NUMPTS (NUMTBL^ 1 )*LEE 

108 NUMT B L-NUMTB L ♦ 1 
GO TO 102 

1101 WRITE ( I B, 1111) LEE 
GO TO 775 

1103 WRITE (IB, 1113 )NUMTBL 
GO TO 775 

1111 FORMAT ( 17H TOO MANY POINTS 18) 

1113 FORMAT (17H TOO MANY TABLES 18) 

69 GO TO (775, 776, 776), NG 

70 IF (MINTBl-MAXm) 71,100,69 

71 DO 73 NAT*MINTBL,MAXTBl 
L4-NUMPTS (NAT )♦ 1 

IF ( ARG3-STG( 14 ) ) 72,74,73 

72 I F ( NAT-MI NTBt )69, 69, 75 

73 CONTINUE 
GO TO 69 

75 15-1 
16*2 
17-14 

101 DO 9/ L8-L5, L6 
L4-NUMPTS (NAT )♦ 1 
L9-L1 (NAT) 

LM-L9+U 
DO 77 LN-1,19 
JJ-UHH 

2626 IF (ARGl-STG(JJ)) 76,78,77 

76 IF (LN-l) 69,69,79 

77 CONTINUE 
GO TO 69 

78 Nl — 1 

GO TO 80 

79 N 1 * ♦ 1 

80 K-L2 (NAT ) 

00 82 1*1, K 
I DATE-LM* I 

IF (ARG2-STG( IDATE ) ) 81,83,82 

81 IF (1-1) 69,69,84 

82 CONTINUE 
GO TO 69 

83 IS— 1 

GO TO 85 

84 IS-*1 

85 I SEQ-LM^12 (NAT)*LN* ( I -1 ) *L9 


02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

0 ? 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 


00530 
00540 
00550 
00560 
00570 
00580 
00590 
00600 
00610 
00620 
00630 
00640 
00650 
00660 
00070 
OC 580 
0069G 
00700 
00710 
00720 
00730 
00740 
00750 
00760 
00770 
00780 
00790 
00800 
00810 
00820 
00830 
00840 
00850 
00860 
00870 
00880 
00890 
00900 
00910 
00920 
00930 
00940 
00950 
00960 
00970 
00980 
00990 
01000 
01010 
01020 
01030 
01040 
01050 
01060 
01070 
01080 
01090 
01100 
oino 
01120 
01130 
01140 
01150 
01160 
01170 
01180 
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J“ I SCQ-L9 
K8»LM*( 1-1) 

K9-L4+LN-1 

IF ( N 1 ♦ I S ) 86,88,91 

86 IF (STG( ISEQ)-999.E20) 87,69,69 

87 FCT«STG( ISEQ) 

GO TO 95 

88 IF (Ml) 89,69,93 

89 IF (AMAXl <STG< ISEQ),STG(J))-999.E20) 90,69,69 

90 FCT"STG( I SEQ)-(STG( I0ATE)-ARG2)*(STG( I SF(l)-STG( J) )/(STG( I DATE) 
l-STG(Ka)) 

GO TO 95 

91 IF (AMAXl <STG(|SEQ),STG(J),STG<ISEQ-1>,STG(J-1))-999.E20> 92, 
1 69,69 

92 FCT«((STG( IDATE)-ARG2>*((STG<JJ)-ARGi)*STGtJ-i>-(STG<K9)-ARGl) 
l*STG(J>>-(STG(K8)-ARG2)«((STG<JJ)-ARGi>*STG(ISEQ-l)-(STG(K9>- 
2ARGi)*STG( ISEQ) ) )/ ( (STG< IOATE)-STG(K8 > >*<STG( JJ)-STG(K9) > > 

GO TO 9b 

93 IF (AMAXl (STG(ISEQ),STG(ISEQ-1))-999.E20) 96,69,69 

94 FCT-STG(ISEQ)-(STG(JJ)-ARG1)*(STG(ISEQ)-STG(ISE(»-1))/(STG(JJ)> 
1STG(K9)) 

95 GO TO (96,98, 99), L8 

96 DUMMY ( 1 )»FCT 

97 HAT-NAT-l . 

98 FCT-DUMMY( 1)-(STG(L7)-ARG3)*<DUMMY(1)-FCT)/(STG<L7)-STG<L4)) 

99 C.0 TO (2000, 3000), NORMAL 
100 MAT-MI NTBL 

74 15-5 
t.6«3 

GO TO 101 

: TNO STINT TABLE LOOK-UP 

U'D 

SUBROUTINE PHI (CG, BS,HKOE,NKDP,NKOFl,NELECT,NPROTT,NFLARE, 
INALPMA.NKDA, PHI TOT) 

I A* 5 
IB-b 

ClASS-f:«/6.0 

SCBACK-BS/5.0 

I’ll I ELI • 0.0 

Pill PRO - 0.0 

PHIFER • 0.0 

PH I ALP • 0.0 

PHI TOT -0.0 

COMPUTE CLECTRON FLUX 

UU 1 I • 1,11 

DUM • FLOAT (I) 

CALL STINT ( OUM , GLASS, 0.0, XKDE , 1 , NGR I PE , NKDE , NKOE ) 

CALL STINT (OUM,SCBACK,0.0,YKOE,1,NGRIPE, NKOE, NKOE) 

CALL STINT <0UM,0.0,0.0,XNE,1,NGR!PE,NEIECT,NELECT) 

PHI ELE - PHIELE ♦ XNE*XKOE ♦ XNE*YKDE 
IT (NGR I PE) 10,1,10 
10 NGR I PC - 5555 
GO TO It 

1 CONTINUE 

C**»* COMPU T E PROTON FLUX 
00 2 I • 1,18 
DUM - FLOAT (I) 

CALL STINT (OUM, GLASS, 0.0,XK0P,l,NGRIPE,NKDP,NKOP) 

CALL STINT (OUM,SCBACK,0.0,YKOP,1,NGRIPE,NKOP,NKOP) 

CALL STINT (OUM,0.0,0.0,XNP,1,NGRIPE,NPROTT,NPROTT) 

PH I PRO - PH I PRO ♦ XNP«XKOP + XNP*YROP 
IF (NGR I PE) 20,2,20 
20 NGRIPE • 6666 
GO TO 11 

2 CONTINUE 



COMPUTE SOLAR FLARE PROTON FLUX 


03 

00350 


DO 3 1 » 1,18 


03 

00360 


DUM - FLOAT (1) 


03 

00370 


CALL STINT (DUM, GLASS, 0 .0, XKDL, 1,NGR 1 PE,NKDF L, NKDFL ) 


03 

00380 


CALL STINT (DUM,SCBACK, 0 , 0, YKDL, 1,NGR 1 PE,NKDFL,NKDFL) 


03 

00390 


CALL STINT ( DUM, 0 . 0, 0 . 0, XNF , 1 , NC.R 1 PE ,NFLARE , NF LARE ) 


03 

00400 


PMIFLR - PHIFLR ♦ XNF *XKDL ♦ XNF*YKl»L 


03 

00410 


IF (NGRIPE) 30,3,30 


03 

00420 

30 

NGRIPE » 7777 


03 

00430 


GO TO 11 


03 

00440 

3 

CONTINUE 


03 

00450 

c***# 

COMPUTE SOLAR FLARE ALPHA PARTICLE FLUX 


03 

00460 


DO L 1 • 1,18 


03 

00470 


DUM • FLOAT (1) 


03 

00480 


CALL STINT (DUM, GLASS, 0. 0, XKDA, 1,NGRI PE, NKDA,NKDA) 


03 

00490 


CALL STINT (DUM, SCBACK,0.0, YKDA, 1,NGRI PE,NKOA,NKOA) 


03 

00500 


CALL STINT (DUM, 0.0,0.0,XNA,1, NGRIPE, NALPHA.NALPMA) 


03 

00510 


PH 1 ALP » PH 1 ALP ♦ XNA*XKDA ♦ XNA«YKDA 


03 

00520 


IF (NGRIPE) UC,u,40 


03 

00530 

40 

NGRIPE • S888 


03 

00540 


GO TO 11 


03 

00550 

4 

CONT 1 NUE 


03 

00560 

C* • • * 

COMPUTE TOTAL FLUX 


03 

00570 


PHITOT • PHIELE ♦ PHIPRO ♦ PHIFI.R ♦ PHIALP 


03 

00580 


WRITE (IB, 100) 


03 

00590 

100 

FORMAT (58H1 SUMMATIONS OF ELECTRON, PROTON, ANO SOLAR 

FLARF 

PARTI C03 

00600 

1LES//89H ELECTRONS PROTONS SOLAR FLARE PROTONS 

S0LA03 

00610 

2R FLARE ALPHA PARTICLES TOTAL//) 


03 

00620 


WRITE ( IB, 101) PHI FLE, PHI PRO, PHIFLR, PHI ALP, PHI TOT 


03 

00630 

101 

FORMAT UX,E12.4,2X,E12.4,Et2.U,15X,E12.k,19X,E12.4) 


03 

00640 


WRITE (IB, 102) 


03 

00650 

102 

FORMAT (54HO NOTE, THESE FLUXES ARE IN EQUIVALENT 1 MF.V ELECTR0NS)03 

00660 


RETURN 


03 

00670 

11 

WRITE ( IP,lQ3)NC»niPF 


03 

00680 


CALL EXIT 


0) 

00690 

103 

FORMAT ( 1 5H0 ERROR IN TLU IS) 


03 

00700 


END 


03 

00710 


SUB ROUT I Nr! DFGRAD ( PH 1 TOT , BOHMS,NCE LIT , MMVTBL ) 


04 

00020 


DIMENSION VVfC< 101 ) , X 1 VEC{ 101 ) 


04 

00030 


COMMON VVFC,XIVFC 


04 

00040 


1 A-5 


04 

00050 


IB-6 


04 

00060 


VVFC(l) • 0*0 


04 

00070 


XlVfCU) • 0*0 


04 

00080 


00 1 1 • 2,101 


04 

00090 


XtVEC(l) • 0*0 


04 

00100 

1 

vvrcci) ■ vvFcn-i) ♦ 0*01 


04 

00110 


DO 3 1 « 1,101 


04 

00120 


CAtL STINT <VVrC<l),0.0,0.0,KIVfCm,l, NGRIPE, NCFLLT, 

NCELLT) 

04 

00130 


IF (NGRIPE) 2,3,2 


04 

00140 

2 

WRITE ( 1 R, 100 ) 


04 

00150 

100 

FORMAT ( 48H0 ERROR IN DEGRADATION ROUTINE TLU OF CELL 

CURVE ) 

04 

00160 


CALL EXIT 


04 

00170 

3 

CONTINUE 


04 

00180 


TLU OF MILLIVOLTS AT ISC VS. TOTAL FLUX AND BASE RtSTIVITY 

04 

00190 


CALL STINT ( PHI TOT, BOHMS.0 .0, Xf'V, 1 ,NOR 1 PE,NMVTBL,NMVTBL ) 

0* 

00200 


IF (NGRIPE) 4,5,4 


04 

00210 

4 

NGRIPE • NMVTBL 


04 

00220 


WRITE ( 1 P, 10 1 )NGRI PE 


04 

00250 

101 

FORMAT ( 25HO ERROR IN TLU, TABLE NO. IS) 


04 

00240 

# 

COMPUTE RR 


04 

00250 

5 

XISCO * XIVEC(l) 


04 

00260 


RR • XMV/XISCO 


04 

00270 

c* * • * 

CHANGE VVEC 


04 

00280 


DO 6 1 * 1,101 


04 

00290 


VVEC(I) - VVEC(I) ♦ XIVEC(I)*RR 


04 

00300 


A-9 


6 

CONTINUE 

04 

00310 

C* * * * 

FIND VOCO 

04 

00320 


00 7 1 - 2/101 

04 

00330 


N* 1 

04 

00340 


IF (XIVEC(D) 8,8/7 

04 

00350 

7 

CONTINUE 

04 

00360 

8 

SLOPE - (XIVEC(N-l)-XIVEC(N-2))/(VVEC(N-l)-VVEC(N-2)) 

04 

00370 


00 9 J - N, 101 

04 

00380 

9 

XIVEC(J) ■ XI VEC(N-1)*(VVEC(J)-VVEC(N-1))*SL0PE 

04 

00390 


VOCO - VVEC(N-l) ♦ XIVEC(N-1)*(VVEC(N)-VVEC(N-1))/(XIVEC(N-1) 

04 

00400 


1 -XIVEC(N)) 

04 

00410 

C * * * • 

COMPUTE Rl FOR TOT* PHI 

04 

00420 


IF (BOHMS - 1.0) 1C, 10,13 

04 

00430 

10 

IF (PHITOT - 1.0E15) 11,12/12 

04 

00440 

11 

IF (PHITOT - 4.0E12) 17,17,18 

04 

00450 

12 

Rl * 3.3 - 0. 167*0. 4342944819*ALOG (PHITOT) 

04 

00460 


GO TO 16 

04 

00470 

13 

IF (PHITOT - 1.0E15) 14,15,15 

04 

00480 

14 

IF (PHITOT - 5.0E12) 17,17,19 

04 

00490 

15 

Rl * 2.806 - 0. 1325*0. 4342944819*ALOG (PH! T OT) 

04 

00500 


GO TO 16 

04 

00510 

17 

Rl - 1.0 

C4 

00520 


GO TO 16 

04 

00530 

18 

Rl » 1.0 - 7. 13/ ((10.0)** 8.31)*((PHITOT-4.0E12)**C.451) 

04 

00540 


GO TO 16 

04 

00550 

19 

Rl = 1.0 - 2.78/((10.0)** 7.34)*((PHITOT-5.0E12)**0.41) 

04 

00560 

c**** 

COMPUTE DELTA 1 

04 

00570 

16 

DELTA 1 * XISCO*(1.0 - Rl) 

04 

00580 

C* * * * 

CHANGE XIVEC 

04 

00590 


DO 20 1 - 1,101 

04 

00600 


X 1 VEC ( 1 ) « X 1 VEC ( 1 ) - DELTA 1 

04 

00610 

20 

CONTINUE 

04 

00620 

c**** 

FIND PZC, POINT OF ZERO CURRENT 

04 

00630 


DO 21 1 - 2,101 

04 

00640 


N « 1 

04 

00650 


IF (XIVEC( 1 ) )22,22,21 

04 

00660 

21 

CONTINUE 

04 

00570 

22 

SLOPE * ( X 1 VEC ( N - 1 ) -X * VEC ( N- 2 ) ) / ( VVEC ( N- 1 ) • VVEC ( N- 2 ) ) 

04 

00680 


DO 23 J « N, 101 

04 

00690 

23 

XIVEC(J) = XIVEC(N-1)*(VVEL(J)-VVEC(N-1))*SL0PE 

04 

00700 


PZC - VVEC(N-1)*XIVEC(N-1)*(VVEC(N)-VVEC(N-1))/'XIVEC(N-1) 

04 

00710 

] 

L -XIVEC(N)) 

04 

00720 

c**** 

COMPUTE RV FOR TOTAL PHI 

04 

00730 


IF (PHITOT - 1.0E14) 25,24,24 

04 

00740 

24 

RV * 1.779 - 0. 0588*0. 4 342 9448 19*ALOG (PHITOT) 

04 

00750 


GO TO 26 

04 

00760 

25 

IF (PHITOT - 3.0E12) 27,27,28 

04 

00770 

27 

RV - 1.0 

04 

00780 


GO TO 26 

04 

00790 

28 

RV *= 1.0 - 0 . 02 2 2 * ( ( 0 . 4 34 29448 19*AL0G ( PH 1 TOT/ 3 . OE 12 ) ) **1 . 67 ) 

04 

00800 

c • * * * 

COMPUTE VOCR 

04 

00810 

26 

VOCR - VOCO*RV 

04 

00820 


DELTAV - PZC - VOCR 

04 

00830 

c**** 

FINAL UPDATE OF VOLTAGES IN ARRAY 

04 

00840 


DO 30 1 - 1,101 

04 

00850 


VVEC(I) - VVcC(l) - DELTAV 

04 

00860 

30 

CONTINUE 

04 

00870 

C* * * * 

COMPUTE ISCR, 1 PMR/VPMR 

04 

00880 

31 

DO 33 1 - 1,101 

04 

00890 


IF ( VVEC (l)> 33,32,32 

04 

00900 

32 

X 1 SCR - X 1 VEC ( 1 ) 

04 

or9io 


GO TO 34 

04 

0(920 

33 

CONTINUE 

04 

00930 

34 

N ■ 1 

04 

00940 


PWR * -1000000.0 

04 

00950 

40 

PWRT - VVEC(N)*XIVEC(N) 

04 

00960 


A- 10 


oooooooooooooo 


41 


44 


42 

43 


48 

45 

102 

103 


105 


106 

46 


F10.6/36H0 
AT MAX PWR 
= F10.6//) 

AMPS//) 


04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

CURREN04 
PT ( VOLT04 
04 


IF ( PWRT-PWR) 42,44,41 
N = N+l 
PWR * PWRT 
GO TO 40 
VPMR * VVEC(N) 

XIPMR = XIVEC(N) 

GO TO 45 

SLOPE = ( X I VEC (N-l ) “X I VEC (N ) )/ ( VVEC(N ) -VVEC (N-l) ) 

VPMR = VVEC(N-l) ♦ O, ,00015 

XIPMR = X I VEC(N-l) - SLOPE* (VPMR - VVEC(N-l)) 

PWRT = V PMR*X I PMR 
IF (PWRT-PWR) 45,45,48 
PWR = PWRT 
GO TO 43 
WRITE (IB, 102) 

FORMAT ( 43H1 SOLAR CELL I -V CURVE IRRADIATION-DEGRADED) 

WRITE ( I B , 1 0 3 ) X I SCR, XJ PMR, VPMR, VO CR 
FORMAT ( 36H0 SHORT-CIRCUIT CURRENT (AMPS) = 

IT AT MAX PWR PT (AMPS) = F10.6/36H0 VOLTAGE 
2S ) = F10.6/36H0 OPEN-CIRCUIT VOLTAGE (VOLTS) 

WRITE (IB, 105) 

FORMAT ( 20H0 DEGRADED !-V CURVE// 17H VOLTS 
DO 46 I * 1,101,2 
WRI TE (IB, 106 ) VVEC ( I ) , X I VEC( I ) 

FORMAT (2F12.6) 

CONTINUE 
RETURN 
END 

SUBROUT I NE STASH (S I G i SC, S I GVOC, D I SC, DV, TNOT, KEY, NGR I PE, NCELLT, 

1VOCO I , ADDT, DELTT) 

DIMENSION WE C ( 1 0 1 ) , X I V E C' ( 10 1 ) ,TEMP( 15 ) , VMAP ( 15 , 101) , X I MAP ( 15 , 1 0 1 ) 05 
D I MENS l-ON VOCT ( 15 ) , X I SCT( 15 ) , VMP( 15 ) ,X I Mp( 15 ) 

DIMENS ION PWRSTG ( 15 ) 

COMMON VVEC, XI VEC 

DECIDE IF INITIALIZATION OR COMPUTATION 
I A = 5 
1 B = 6 

NGR I PE=0 

I F( KEY) 2000, 3000, 1000 

INITIALIZATION - OBTAIN I AND V INPUT VECTORS 
************ ********* ******* * ** * ************** ***************** 

IN FIRST ENTRY 
SIGISC IS DISC 
SI GVOC IS THETA 
DISC IS AVPMO 
DV I S A I PMO 
TNOT IS TNOT 
KEY IS -1 
NGR I PE IS NGR I PE 
NCELLT IS NCELLT 
VOCO? i $ VOCOI 
ADPT IS ADDT 
DELTT IS DELTT 

2000 CALL SLITET (1, K001FX) 

GO TO (2005,5000), K001FX 

5000 VVEC ( 1 ) =0 . 0 
X I VEC ( 1 ) = 0 . 0 
DO 2001 1*2,101 
X I VEC ( I )*0.0 

2001 VVEC ( I ) “ VVEC(l-l) ♦ 0.010 
DO 2002 I - 1,100 

CALL STINT (VVEC( I ) , 0 . 0 , 0 . 0, X I VEC ( I ) , 1, NGR I PE, NCELLT, NCELLT) 

I F( NGR I PE) 2002, 2002, 2003 

2002 CONTINUE 
GO TO 2005 

2003 I F ( X I VEC ( I -1 ) ) 2005 , 2005, 12003 


04 

04 

04 

04 

04 

04 

04 

04 

05 
05 


05 
05 
05 
05 
05 
0 5 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 
05 


A-ll 


00970 
00980 
00990 
01000 
01010 
01020 
01030 
01040 
01050 
01060 
01070 
01080 
01090 
01100 
OHIO 
01120 
01130 
01140 
01150 
01160 
01170 
01180 
01190 
01200 
01210 
01220 
01230 
01240 
00020 
000 30 
00040 
00050 
00060 
00070 
00080 
00090 
00100 
00110 
00120 
00130 
00140 
00150 
00160 
00170 
00180 
00190 
00200 
00210 
00220 
00230 
00240 
00250 
00260 
00270 
00280 
P0290 
00300 
00310 
00320 
00330 
00340 
00350 
00360 
00370 
00380 
00390 


12003 

RETURN 

05 

oouoo 

C 

EXTRAPOLATE INPUT CURVE 

05 

00410 

C 

LOCATE FIRST ZERO ELEMENT IN XIVEC 

C5 

00420 

2005 

DO 2107 1 - 2,101 

05 

00430 


N»l 

05 

00440 


1 F ( X 1 VEC( 1)12106,2106,2107 

05 

00450 

2107 

CONTINUE 

05 

00460 

2106 

SLOPE«(XIVEC(N-l)-XIVEC<N-2))/(WEC(N-l)-VVEC(N-2)) 

05 

00470 

C 

LOCATE POINT OF ZERO CURRENT FOR '•NOERRADED CURVE 

05 

0C480 


PZCU • VVEC(I) - X 1 V EC ( 1 l/SLOPE 

05 

00490 


DO 2108 J » N , 1 0 1 

05 

00500 

2108 

XIVEC(J)=XIVEC(N-l)*(VVEC(J)«VVEC(N-l))*SLOPE 

05 

00510 

C 

00 CURRENT DEGRADATION WITH GAMMA, TNOT, THETA 

05 

00520 


XISC « XIVEC(l) 

05 

00530 


GAMMA=S (RISC 

05 

00540 


DELTAI*(1, 0 -GAMMA ) »X 1 SC 4 

05 

00550 


REC 1 P*1 , 0/ GAMMA 

05 

00560 


THETA=S 1 RVOC 

05 

00570 



PFLTAV=(TNOT+273.16 ) *0 . 8614E-04* ALOG 

< REC 1 ? ) 

05 

00580 


OELTAV * nr.LTAV ♦ <1.0 - THETA)*VOCOI 



05 

00590 


DO 2600 1 » 1,101 



05 

00600 


XIVLCd ) • XIVFC(I) - DELTAI 



05 

00610 


IF (XIVEC(D) 2703,2705,2600 



05 

0C620 

2705 

MX » 1 - 1 



05 

00630 


IOU « »♦ 1 



05 

UP 64 0 


00 2 G2 2 t * LOO # 101 



05 

OQGSO 

2622 

XIVEC(L) « Xlvrc(L) - OEtTAI 



05 

C0b( 0 


<;o io ? no*. 



05 

096 70 

2600 

LUNt INtll 



05 

00660 

2006 

CONT 1 MO? 



05 

00690 

C 

LUCA If: I’0 1 N 1 01 ZFRO CURRENT FOR DEGRADED 

CURVE 

05 

00700 


S l.ui'L ’ (XlVff ( K 1 X ) - XIVEC(KIX>l))/(VVECOCIX) - VVEC(KIX-l)) 

05 

00710 


PZU* • VVfUMX) - XIVCGUIXl/fl.OPF 



05 

00720 


UEITAV • (IVCU - PZCD) • OCITAV 



0* 

00 7 30 


t'fi 2 701. 1 • 2,101 



05 

00740 

2704 

vvrcd) * vvrrm ♦ dfltav 



05 

00750 

C 

C'lltniftT I>* M:Ar A 1 ION IS NON COMPLETE 



C5 

0076C 

c 

r »>^ ON !*• n r ff irs PE S 1 STANCE DEGRADATION 


r r 

0 r 7?ft 

c 

STOP* MW: AVPMO AND FETCH NCK BATCH 

OF VARIARLFS 

95 

D0780 


Air?'0«HV 



05 

CC 790 


AVPMG*0 1 C C 



05 

00800 


RETURN 



05 

00810 

£***•< 


t • t • < 


C5 

OCC20 

c 

IN SrCIHID f At 1 



05 

00830 

c 

r V\\ «C 1°) ( V 



05 

00849 

c 

SIOVOC IS DISC 



05 

00850 

c 

List i$ staisr 



05 

00860 

c 

pv is s invar 



05 

00870 

c 

tnot is t ni r 



05 

00880 

c 

KEY IS 0 



05 

00800 

c 

NGRtPE IS NVtlPE 



05 

00900 

c 

NCCLLT IS NCflLT 



05 

00910 

c 

VOCOI IS OEHFAC 



OS 

00920 

c 

A DOT IS AIN'T 



05 

00930 

c 

DftTT IS DfLTT 



05 

00940 

3000 

ALFA*S IS ISC 



05 

00950 


DENFAC«VCCOI 



C5 

00960 


TM«tnsc 



05 

0C970 


TKV«DV 



05 

00980 


RES«((l.0-ALFA)»AVPMO/AIPMO) 



05 

00990 


DO 3001 ! • 1,10* 



05 

01000 

3001 

VVEC ( 1 )«VVf Cd ) *X 1 VEC ( 1 )«RFS 



05 

01010 

c 

SERIES RESISTANCE DEGRADATION IS NOW 

COMPUTE 

05 

01020 


A-12 



c 

DETERMINE VOC - LOCATE FIRST NEGATIVE ELEMENT IN XIVEC 

05 

01030 


DO 3007 I • 2,101 


05 

01040 


N-l 


05 

01050 


1 F (X 1 VEC ( 1 >>3006,3006,3007 


05 

01060 

3007 

CONTINUE 


05 

01070 

3006 

VOC-VVEC(N-1)+XIVEC(N-1)*(VVEC<N)-VVEC(N-1>)/(XIVEC(N-1)-XIVEC(N))05 

01080 

C 

FILL TEMPERATURE VECTOR 


05 

01090 


IF (DELTT) 3012,3012,3013 


05 

01100 

3012 

ADDT * 120.0 


05 

OHIO 


DELTT » 20.0 


05 

01120 

3013 

CONTINUE 


05 

01130 


TEMP(l) - TNOT «■ ADDT 


05 

01140 


V0CT(1)» VOC-ADDT*TKV 


05 

01150 


DO 3202 1-2,15 


05 

01160 


TEMPO ) • TEMPO-1) - DELTT 


05 

01170 

3202 

VOCT ( 1 )=VOC-(TEMP( 1 )-TN0T)*TKV 


05 

01180 

C 

EXPAND 1 -V CURVE INTO FAMILY OF CURVES FOR 

15 TEMPERATURES 

05 

01190 

C 

OUTER DO FOR EACH TEMPERATURE ON INDEX 1 


05 

01200 


DO 3040 1-1,15 


05 

01210 

C 

FIRST INNER DO FOR VOLTAGES ON INDEX J 


05 

01220 


T2T1-TEMP ( 1 )-TNOT 


05 

01230 


DELTAV-TKV-T2T1 


05 

01240 


DO 3204 J - 1,101 


05 

01250 

3204 

VMAP (l,J)-VVEC(J)-DELTAV 


05 

01260 


XK 1 T-TK 1 *T2T1*S 1 GVOC 


05 

01270 


DENOM-VOC-DENFAC-TKV-T2T1 


05 

01280 

C 

SECOND INNER DO FOR CURRENTS ON INDEX J 


05 

01290 


DO 3205 J - 1,101 


05 

01300 


Z-VMAPO, J)/DENOM 


05 

01310 


Z6 » Z**8 


05 

01320 


UMZ6-1 . 0-Z6 


05 

01330 


IF(UMZ6)33Q0, 3205, 3205 


05 

01340 

3300 

UMZ6-0.0 


05 

01550 

3205 

XIMAP( l,J)=XIVEC(J)+XKIT*UMZ6 


05 

01360 

3040 

CONTINUE 


05 

01370 

C 

FIND V AND 1 AT MAXIMUM POWER POINT FOR EACH TEMPERATURE 

05 

01380 


DO 4005 1-1,15 


05 

01390 


N-4 


. 05 

01400 


PWR- -1000000.0 


05 

01410 

4000 

PWRT-VMAP ( 1 , N ) *X 1 MAP ( 1 , N ) 


05 

01420 


1 F ( PWRT-PWR) 4002,4004,4001 


05 

01430 

4001 

N-N+l 


05 

01440 


PWR-PWRT 


05 

01450 


GO TO 4000 


05 

01460 

4004 

VMP( 1 )*VMAP( l,N) 


05 

01470 


X 1 MP( 1 )=XIMAP( 1 ,N) 


05 

01480 


GO TO 4005 


05 

01490 

4002 

SLOPE-C X 1 MAP ( 1 , N-l) -X 1 MAP( 1 ,N )>/(VMAP( 1 , N >- 

VMAP( 1 ,N-1) ) 

05 

01500 

4003 

VMP( l)-VMAP( l,N-l)+0. 00015 


05 

01510 


X 1 MP (1 ) —XI MAP C 1 , N-l ) -SLOPE* (VMP( 1 >-VMAP( l,N 

-1)> 

05 

01520 


PWRT=VMP( 1 > * X 1 M P ( 1 > 


05 

01530 


IF (PWRT-PWR) 4005,4005,4008 


05 

01540 

4008 

PWR-PWRT 


05 

01550 


GO TO 4003 


05 

01560 

4005 

CONTINUE 


05 

01570 


DO 100 1 - 1,15 


05 

01580 


PWRSTG ( 1 ) « VMP( 1 )*XIMP< 1 ) 


05 

01590 

10Q 

CONTINUE 


05 

01600 

4006 

DO 4009 1 - 1,15 


05 

oi6 10 


DO 4010 N - 1,101 


05 

01620 


IF ( VMAP ( 1 ,N) >40 10,4007,4007 


05 

01630 

4007 

XISCT(I) - X IMAP( 1 ,N) 


05 

01640 


GO TO 4009 


05 

01650 

4010 

CONTINUE 


03 

01660 

4009 

CONTINUE 


05 

01670 

C 

EDIT RESULTS 


05 

01680 


A-13 



WRITE (16,7000) 

WRITE ( IB,7001)(TEMP(J),J-1,15) 

WRITE (IB,7002)(XIMP(J),J-1,15) 

WRITE (IB,7003)(VMP(J),J-i,15) 

WRITE (IB, 700 6 ) ( VOCT ( I ) , I *1, 15 ) 

WRITE (IB,7007)(XISCT(I ), 1*1,15) 
WRITE (»B,7008)( PWRSTG( I), 1*1,15) 

00 2040 d-1,30,2 

WRITE (18,7004)*), (X IMAP ( I, J), 1*1,15) 

2040 WRITE ( I B, 7005 ) J, (VMAP( I ,*)) , 1*1,15) 
WRITE (16,7000) 

WRITE (IB,7001)(TEMP(*J),J-1,1S) 

WRITE (IB,7002)(X)MP(J),J-1,15) 

WRIT* ( IB,7003)(VMP(J), J-1,15) 

WRITE (16,7006) (VOCT (I), 1*1, 15) 

WRITE (IB, 7007)(XISCT(I), 1-1,15) 
WRITE ( I B, 7008 ) (PWRSTG( I), 1-1,15) 

00 2041 J - 31,60,2 

WRITE ( I B, 7004 )J, (X IMAP( I, J), 1-1,15) 

2041 WRITE (IB,7005)J, (VMAP( l,J), 1-1,15) 
WRITE (IB, 7000) 

WRITE < IB,7001)(TEHP(J),J»1,15> 

WRITE (IB,7002)(XIMP(J),J-1,15) 

WRITE (16,7003) ( VMP(J ) , *1-1, 15 ) 

WRITE { IB,7006)(VOCT( I ), 1-1, 15) 

WRITE (IB,7007)(XISCT( I ), 1-1,15) 
WRITE ( IB, 7008) (PWRSTG( I), 1-1,15) 

00 2042 J - 61,90,2 

WRITE (18,7004)*), (X«MAP( I, J), 1-1,15) 

2042 WRITE ( I B, 7005 ) J, (VMAP( l,J), 1-1, IS) 
WRITE (IB, 7000) 

WRITE (18, 70P1)(TEMP(*J), *J-1,15) 

WRITE ( I P.,7002 )(XIMP(J),J-1, 15) 

WRITE (IB, 7003) ( VMP( J ) , J-l , IS ) 

WRITE ( I B, 7006 ) (VOCT ( I ) , I " 1,15) 
WRITE ( I B, 7007 ) (X I SCT( 0,1-1,15) 
WRITE ( I B, 7008 ) (PWR$TG( I), 1-1,15) 

00 2043 J - 91,99,2 

WRITE ( 16,7004 )J, (XIMAP( I ,J>, I -1,15) 

2043 WRITE ( I P, 7005 ) J, (VMAP( I , J), I -1, 1$ > 
RETURN 


05 01690 
05 01700 
05 01710 
05 01720 
05 01730 
05 01740 
05 01750 
05 01760 
05 01770 
05 01780 
05 01790 
05 01800 
05 01810 
05 61820 
05 01830 
05 01840 
05 01250 
05 01860 
05 01870 
05 01880 
05 01890 
05 01900 
05 01910 
05 01920 
05 01930 
05 01940 
05 01950 
OS 01960 
05 01970 
C5 01980 
05 01990 
05 02000 
05 02010 
05 02020 
05 02030 
05 02040 
05 02050 
05 02060 
05 02070 
05 02080 
05 02090 


EDIT FORMATS 
FORMAT ( 1041*1 


EOIT OF TEMPERATURE AND DEGRADATION CORRECTED SOLAKOS 


l CFLl l-V CURVES FOR SOLAR ARRAY SYNTHESIS PROGRAM) 
FORMAT ( 14H0 TEMPERATURES15F7.0) 

FORMAT U4HQ I -MAX PWR 15F7.4) 

FORMAT (1411 V-MAX PWR 15F7.4) 

FORMAT ( 4HO 1(12,810 15F7.4) 

FORMAT (4H V(I2,8H) 15F7.4) 

FORMAT (14,1 VOLTS OC 15F7.4) 

FORMAT (14H AMPS SC 1SF7.4) 

FORMAT (14H MAX POWER 15F7.4) 


LOOK UP CURRENT AND MAX PWR POINT GIVEN VCELL ♦ TEMPERATURE 05 
VCELL-DI SC 05 
HELHOT-OV 05 
NK - (TEMPI 1 ) - HELHOT)/DELTT 05 
NHI-NK-1 05 
NLO-NK*2 05 
IF (KEY - 989) 1334,1333,1334 05 
FOR KEY - 989 INTO SIGISC PUT OPEN CIRCUIT VOLTAGE AT TEMPETATURE 05 
SIGISC - VOCT(NLO) ♦( (OV-TEMP(NLO) )/0ELTT)*( VOCT(NHI >-VOCT(NLO) ) 05 
RETURN 05 
IF (KEY - 988) 1335,1336,1335 05 
FOR KEY - 988 INTO SIGISC PUT VOLTAGE AT THE GIVEN CURRENT 05 
00 1402 J > 2,101 05 
N - *1 05 


A- 14 





-~-T» 




X IMAP(NLO, J) ) 1402, '402,1403 


(VMAP(NL0,N> - 


IF (VCEll 

1402 CONTINUE 

1403 XILO - VMAP(NLO,N-l) ♦ (VCELL - XIMAP(NLO,N-l) >« 

1VMAP(NL0,N-1) )/ (X IMAP(NLO,N)-XIMAP(NLO,N-1) ) 

DO 1422 J ■ 2,101 

N « J 

IF (VCELL- XIMAP(NHI,J))1422, 1422, 1423 

1422 CONTINUE 

1423 CONTINUE 

X IH I » 7MAP(NHI,H-1) ♦ (VCELL - X IMAP(NH I, N-l ) )* (VMAP(NHI,N) 
1VMAP(NHI,N-1))/(XIMAP(NHI,N) - X IMAP(NHI ,N-1 ) ) 

GO TO 1210 
1335 CONTINUE 

I F (VCELL-VMAP(NLQ, 1 ) >1100, 1100,1001 
1100 XILO*XIMAP(NLO,l) 

GO TO 1010 

1001 00 1002 J « 2,101 
N«J 

IF(VCELL-VMAP(NLO,J>) 1003, 1003, 1002 

1002 CONTINUE 

1003 X I LG*X IMAP(NLO,N-l )♦ ( VCELL*VMAP(NLO,N-l ) )* (X IMAPtNLO, N)-X IMAPCNLO, 05 

IN-1) )/ ( VMAP(NLO,N)-VMAP(NLO,N~l) ) 05 

1010 I F ( VCELL- VMAP(NH I , 1 )) 1200, 1200, 1020 05 

1200 X IM I “X IMAP(NH 1,1) 05 

GO TO 1210 05 

1020 DO 1222 J - 2,101 05 

N»*J 05 

I F(VCELL-VMAP(NH I, J) >1023,1023, 1222 05 

1222 CONTINUE 05 


05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 

05 


1023 XIHI-XIMAP(NHI,N-1)4(VCELL-VMAP(NHI,N-1>>*(X*MAP(NHI,N)-XIMAP(NHI,05 


IN-1 ) )/ ( VMAPCNH 1 ,N)-VMAP(NHI , N*1 ) ) 05 

1210 FAKTOR»(HELHOT-TEMP(NLO))/<TEMP(NHI)-TEMP(NLO)> 05 

C X I CELL PUT INTO SIGISC 05 

$ IGISC"X I LO+(X IHI -X I LO)*FAKTOR 05 

RETURN 05 

END 05 


02350 

02360 

02370 

02380 

02390 

02400 

02410 

02420 

02430 

02440 

02450 

02460 

02470 

02480 

02490 

02500 

02510 

02520 

02530 

02540 

02550 

02560 

02570 

02580 

02590 

02600 

02610 

02620 

02630 

02640 

02650 

02660 

02670 

02680 

02650 

02700 


A-15 


1 


mh- 

■ 


Is 


Appendix B 


071768 o 
SAMPLE TABLES 


1 OHM-CM %N IMBUS D< 


00 


05- 


0.0 

0.25 

0.30 

0.35 

0.37 

0.39 

0.40 

0.41 

0.42 

01 

0.1420 

0.1420 

0.1412 

0.1412 

0.1409 

0.1400 

0.1397 

0.1390 

0.1380 

02 

0.43 

0.44 

0.45 

0.46 

0.47 

0.48 

0.49 

0.50 

0.51 

03 

0.1369 

0.1355 

0.1332 

0.1312 

0.1289 

0.1256 

0.1210 

0.1161 

0.1088 

04 

0.52 

0.53 

0.54 

0.55 

0.56 

0.57 

0.590 

0.600 

1.0 

05 

0.1005 

0.0925 

0.0825 

0.0700 

0.0550 

0.0388 

0.0 

-0.0125-0.012506 

6800023601 RELATIVE 

ISC V$ 

INCIDENCE ANG'.ESKELLY COSINE< 

00 

0.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40,0 

01 

1.0 

0.99 

0.98 

0.965 

0.94 

0.905 

0.86 

0.81 

0.75 

02 

4S.0 

50.0 

60.0 

65.0 

70.0 

75.0 

80.0 

85.0 

90.0 

03 

0.69 

0.63 

0.48 

0.385 

0.30 

0.20 

0.12 

0.04 

0.0 

04 

270.0 

275.0 

280.0 

285.0 

290.6 

295.0 

300.0 

310.0 

315.0 

05 

0.0 

0.04 

0.12 

0.20 

0.30 

0.385 

0.48 

0.63 

0.69 

06 

320.0 

325.0 

330.0 

335.0 

340.0 

345.0 

350.0 

355.0 

360.0 

07 

0.7S 

0.81 

0.86 

0.905 

0.94 

0.965 

0.98 

0.99 

1,0 

08 

6800031 

101 NIMBUS ORBITAL 

ELECTRONS 1 YR 

600NM- 

1970 


00 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

01 


2.924121.104125.484117.3 41037. 24107. 814101. 614105. 40 491.72 4902 


10*0 11*0 
5.91 480.0 

07-50-6800041801 NIMBUS ORBITAL PROTONS 1 YR 600NM-1970 

1.0 2.0 5.0 4.0 5.0 6.0 7.0 8.0 9.0 

2. 194102. 194101. 4 6410 1.82ft 101 .284)01. 1 4101.1 4109.12 497.3 

10.0 11.0 12.0 13.0 14.0 15.0 16.0 17.0 18.0 

5.48 495.48 494.74 493.28 492.56 491.82 4912.8 696. 9U 690.0 

NIMBUS SOLAR FLARE PROTONS 1YR 600NM 

3.0 4 v 0 5.0 6.0 7.0 

0.0 0.0 0.0 0.0 0.0 

12.0 13.0 14.0 15.0 16.0 

0.0 0.0 0.0 0.0 0.0 


05-13-6800051801 

1.0 2.0 

0.0 0.0 

10.0 11.0 

0.0 ' 0.0 

05-13-6800061801 

1.0 2.0 

0.0 0.0 

10.0 11.0 

0.0 0.0 

09-20-6800071107 


NIMBUS SOLAR FLARE ALPHA PARTICLES 


8.0 
0.0 
17.0 
0.0 
1YR 600NM 


9.0 

0.0 

18.0 

0.0 


3.0 

0.0 

12.0 

0.0 

OAMAGE 

3.0 

0.18 

0.13 

0.08 

0.03 

0.0 

0.0 

0.0 


1.0 2.0 

0.0 0.01 0.06 

1.0 0.0 0.03 

1.5 0.0 0.02 

3.0 0.0 0.0 

6.0 0.0 0.0 

9.0 0.0 0.0 

200.0 0.0 0.0 

10.0 11.0 

0.0 7.4 7.8 

1.0 7.3 7.8 

1.5 7.3 7.8 

3.0 6.8 7.7 

6.0 6.5 7.5 

9.0 5.85 7.0 

200.0 0.0 0.0 

09-20-6800081808 DAMAGE 

1.0 2.0 3.0 


4.0 

0.0 

13.0 

0.0 

FACTORS 

4.0 
0.38 
0.3 
0.2 
0.1 
0.02 
0.0 
0.0 


5.0 6.0 

0.0 0.0 

14.0 15.0 

0.0 0.0 

FOR ELECTRONS 


5.0 

1.3 

1.15 

1.02 

0.75 

0.47 

0.25 

0.0 


6.0 

2.9 

2.7 

2.5 

2.05 
1.55 
1.1 

, 0.0 


7.0 

0.0 

16.0 

0.0 

7.0 

4.35 

4.15 
3 92 
3. >8 
2.85 

2.15 

0.0 


8.0 

0.0 

17.0 

0.0 

8.0 

5.5 
5.3 
5.15 

4.6 
4.1 
3.5 
0.0 


9.0 
O.C 

18.0 

0.0 

9.0 

6.5 
6.15 

6.1 
5.7 
5 3 

4.6 

0.0 


0.0 

0.5 

1.0 

1.5 

3.0 

6.0 

9.0 

200.0 

0.0 


20000.010000.07000.0 


FACTORS 
4.0 


0.0 

0.0 

0.0 

0.0 

0.0 

3.0 

0.0 

10.0 


0.0 

0.0 

0.0 

0.0 

Q.O 

0.0 

0.0 

11.0 


6000.0 
0.0 
0*0 
0.0 

0.0 
0.0 
0.0 
12.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

13.0 


FOR PROTONS 
5.0 


0.0 

0.0 

0,0 

0.0 

0.0 

14.0 


03 

04 
00 
01 

6902 

03 

04 
00 
01 
02 

03 

04 
00 
01 
02 

03 

04 
00 
01 
02 
03 

Ot* 
0- 
06 

07 

08 

09 

10 
11 
12 

13 

14 

15 

16 
00 


6.0 

7.0 

8.0 

9.0 

01 

3900.0 

3500.0 

3100.0 

2800.0 

02 

3900.0 

3500.0 

3100.0 

2800.0 

03 

3700.0 

3500.0 

3100.0 

2800.0 

04 

2000.0 

3100.0 

3100.0 

2800.0 

05 

0.0 

0.0 

200.0 

1400.0 

06 

0.0 

0.0 

0.0 

0.0 

07 

0.0 

0.0 

0.0 

0.0 

08 

0.0 

0.0 

0.0 

0.0 

09 

15.0 

16.0 

17.0 

18.0 

10 

2500.0 

2500.0 

2300.0 

1000.0 

11 


B-l 


> 


: rf J* 




t/ 






at 


XT' 


0,5 

2700.0 

2600.0 

2500.0 

2500.0 

2500.0 

2500. 0 

2500.0 

2300.0 

1000.0 

12 

1.0 

2700.0 

2600.0 

2500.0 

2500.0 

2500.0 

2500.0 

2500.0 

2300.0 

1000.0 

13 

1.5 

2700,0 

2600.0 

2500.0 

2500,0 

2500.0 

2500.0 

2500.0 

2300.0 

1000.0 

14 

3.0 

2000.0 

2100.0 

2100.0 

2100.0 

2000.0 

2000.0 

2000,0 

2000.0 

1000.0 

15 

6.0 

0.0 

0.0 

100.0 

1000.0 

1400.0 

1500.0 

1800,0 

2000.0 

1000.0 

16 

9.0 

0.0 

0.0 

0.0 

0,0 

0.0 

110.0 

1S0Q.0 

2000.0 

1000.0 

17 

200.0 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

18 

09-23- 

6800091807 DAMAGE FACTORS 

FOR ALPHA PARTICLES 



00 


1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

01 

0.0 

20000.1 

015000.1 

014500. 1 

013700. 

012500.1 

011500.1 

011000.1 

010400, 1 

09800.0 

02 

1.0 

lOOOO.i 

013000,1 

014000.) 

013700. 

012500.' 

011500.1 

011000.010400.1 

09800.0 

03 

1.5 

0.0 

0.0 

7000.0 

11000. 

012000. 

011500.1 

011000.1 

010400.1 

09800.0 

04 

3.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2500.0 

5200.0 

7200.0 

7900,0 

05 

6.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

06 

9.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

07 

200.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

08 


10.0 

11.0 

12.0 

13.0 

14.0 

15.0 

16.0 

17.0 

18,0 

09 

0.0 

9500.0 

9500.0 

9500.0 

9500.0 

9400.0 

9200.0 

8600.0 

7000.0 

4000.0 

10 

1.0 

9500.0 

9500.0 

9500.0 

9500.0 

9400,0 

9200.0 

8600,0 

7000.0 

4000.0 

11 

1.5 

9500.0 

9500.0 

9500.0 

9500.0 

9400.0 

9200.0 

8600.0 

7000.0 

4000.0 

12 

3.0 

8000.0 

7800.0 

7700.0 

7600.0 

7600.0 

7600.0 

7700.0 

7000.0 

4000.0 

13 

6.0 

1700.0 

3650.0 

5300.0 

5700.0 

6300.0 

7000.0 

7400.0 

7000.0 

4000.0 

14 

9.0 

0.0 

0.0 

0.0 

1600.0 

4000.0 

6400.0 

7100.0 

7000.0 

4000.0 

15 

200.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

16 

09-20- 

6600102302 NMVTBL CURVE SHAPE VOLTS VS PHI 



00 


1.0 &132.5 4136.0 4131.0 4142.0 4143.0 4144.0 4145.0 41401 


1.0 0.0 0.0 0.001 0.0032 0.005 0.0081 0.0103 0.012 0.0134 02 

10.0 0.0 0.0 0.0 0.0012 0.0024 0.004 0.0051 0.0058 0.0064 03 

6.0 4148.0 4141.0 4151.5 4152.0 4152.5 4153.0 4153.5 4154.0 41504 

1.0 0.0145 0.0163 0,017? 0.0201 0.0219 0.0232 0.0244 0.0253 0.0261 05 

10.0 0.0068 0.007? 0.0085 0.0104 0.0118 0.0130 0,014 0.0149 0.0156 06 

5.0 4156.0 4158.0 4151.0 4161.0 417 07 

1.0 0.0275 0.0286 0.0304 0.0317 0.0457 08 

10.0 0.0168 0.0129 0.0195 0.0207 0.0329 09 

07*17-6800110201 ZERO INCIDENCE ANGLE 00 

0.0 120.0 01 


NIMBUS SOLAR ARRAY TEMP VS TIME PROfilLE 


0.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 

01 

•60.0 

-31. 

0 -11,0 

5.0 

17.0 

27.0 

36.0 

42.0 

47.0 

02 

45.0 

50.0 

55.0 

60.0 

65.0 

70.0 

75.0 



03 

51.0 

52.0 

51.0 

48.0 

46.0 

44.0 

42,0 



04 

.1 3 

MONTH 

FLUX IS 

7.9EXP13 








1 20.0 
2 5.0 

3 0.00007 

4 0.0022 

5 100.0 

6 100.0 

7 100.0 

8 95.0 

9 99.0 

10 100.0 

11 0.47 

12 0.1289 

13 0.590 

14 1O0.O 

15 28.0 

16 1.8 

17 10.0 

18 32.0 

19 79000600000000.0 

20 6.0 

21 15.0 

22 1.0 

23 1.0 






Amt# IT S' 


m Ji i / i 


92.0 

18.0 


11.0 

11.0 


12.0 

12.0 


24 4.0 
99 
2 

4.0 

2.0 

RUN NO. 2 6 MONTH FLUX IS 1.58EXP14 

5 103.0 

6 101.0 

7 103.5 

14 101.0 

19 158000000000000,0 
24 3.0 



B 







